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' Large- multi- matrix loop equations are formulated as quadratic difference equa- 

tions in concatenation of gluon correlations. Though non-linear, they involve highest 
rank correlations linearly. They are underdetermined in many cases. Additional lin- 
ear equations for gluon correlations, associated to symmetries of action and measure 
Q^' are found. Loop equations aren't differential equations as they involve left annihila- 

O . tion, which doesn't satisfy the Leibnitz rule with concatenation. But left annihilation 

is a derivation of the commutative shuffle product. Moreover shuffle and concate- 

■ nation combine to define a bialgebra. Motivated by deformation quantization, we 
^ . expand concatenation around shuffle in powers of q, whose physical value is 1. At 

^ I zeroth order the loop equations become quadratic PDEs in the shuffle algebra. If 

- - ' the variation of the action is linear in iterated commutators of left annihilations, 

these quadratic PDEs linearize by passage to shuffle reciprocal of correlations. Re- 
markably, this is true for regularized versions of the Yang-Mills, Chern-Simons and 
Gaussian actions. But the linear equations are underdetermined just as the loop 
equations were. For any particular solution, the shuffle reciprocal is explicitly in- 
verted to get the zeroth order gluon correlations. To go beyond zeroth order, we find 
a Poisson bracket on the shuffle algebra and associative (7 -products interpolating 
between shuffle and concatenation. This method, and a complementary one of de- 
forming annihilation rather than product are shown to give over and underestimates 
for correlations of a gaussian matrix model. 
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1 Introduction 



1.1 General Remarks 

Approximation methods in physics are often usefully organized as an expansion in a dimension- 
less parameter. As is well known, at first sight, quantum Yang-Mills theory does not have any 
such expansion parameter since the dimensionless coupling of the classical theory is deter- 
mined in terms of the ratio ^ where is the momentum transferred to a hadronic system 
by an external (say electroweak) current. A (say Aqcd) is the dimensional parameter arising 
via dimensional transmutation and renormalization. The success of an expansion in inverse 
(logarithmic) powers of ^ is, however, crucially dependent on the asymptotic freedom of the 
theory for large values of this parameterpp. Thus, this expansion (perturbative QCD), which is 
the analogue of the Born approximation of atomic physics, though spectacularly successful at 
high momentum transfers, is not particularly useful to describe 'intrinsic' properties of hadrons 
in the absence of an external probe transferring a large momentum Q^. 

What about Ti as an expansion parameter for quantum Yang-Mills theory around its classical 
limit? This is a bad starting point, since all variables, not just gauge-invariant ones, stop 
fluctuating in this limit. Since fi can be absorbed into , the 'loop' expansion in powers of 
h around the trivial solution to classical Yang-Mills theory is the same as perturbative QCD. 
Thus, it is useful only at high momentum transfers. 

As observed by 't Hooft[SI, 1/-^ the gauge group SU{N) is an expansion parameter 
for quantum Yang-Mills theory, holding A = g'^N fixed. There are many indicationsj^ that 
N ^ oo is a good approximation to the quantum theory. Moreover, it is a classical limit where 
fluctuations in gauge-invariant variables alone vanish. Despite effort, the 1/N expansion has 
not been as quantitatively successful as perturbative QCD was in the high energy regime. The 
success of the loop expansion lay in the availability of explicit solutions to classical Yang-Mills 
theory around which to expand (eg. flat connections, Euclidean instantons). By contrast, we 
don't know the zeroth order solution of large N Yang-Mills theory around which to perform 
a 1/N expansion. Difficulties are encountered in each of the many ways of formulating the 
large limit of Yang- Mills theory: summing an infinite class of planar diagrams [HI, solving 
the Makeenko-Migdal equations for Wilson loops [31 [HI [Tj or solving the factorized Schwinger- 
Dyson equations for gluon correlations. It would really help to have yet another dimensionless 
expansion parameter, to organize an approximate solution of A = oo Yang- Mills theory. 

The strategy of looking for an expansion parameter over and above 1/A has found success 
in maximally super-symmetric Yang-Mills theory. In some sectors of the AA = 4 theory, an 
expansion around small values of the ratio of 't Hooft coupling to square of ii-charge ( -p ) 
has been developed 0. An analog of this for the non-supersymmetric theory would be useful. 
But since there is no such obvious expansion parameter, we will invent one based on deeper 
mathematical structures of the theory. 

Inspiration for a possible approximation comes from atomic physics, as emphasized by 
Rajeev[n|- The Hartree-Fock approximation for many-electron atoms is analogous to the N ^ oo 
limit of Yang- Mills theory, since it can be formulated as the limit in which the number of replicas 
of each electron (A^) tends to infinity . In general, the Hartree-Fock equations are difficult 
to solve since they involve the electron density matrix, which is a projection operator. However, 
after A — > oo it is possible to take a semiclassical limit based on deformation quantization. 
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These limits do not commute. At zeroth order this leads to the Thomas-Fermi non-linear ODE 
whose solution gives a good first approximation to the charge density of a many-electron atom 
[0]. Can something similar work for large N Yang-Mills theory? 

The approximation method studied in this paper is based on the observation that even in 
the 'classical' large- limit, the equations of matrix models and Yang-Mills theory still involve 
non-commutative concatenation products. It should be possible to take a further 'classical' 
limit, where they are approximated by commutative products by analogy with deformation 
quantization. In our case, the parameter controlling this further classical limit is a deformation 
parameter whose physical value is q = 1. 

Another lesson from the formulation of Hartree-Fock theory as the limit of a large number 
of electron replicas, is that the physical value of an expansion parameter need not be small 
for the expansion to be practically successful. Indeed, the physical number of replicas of the 
electron is A^ = 1 and yet, Hartree-Fock, which corresponds to N = oo, provides a good first 
approximation as part of a 1/A^ expansion! Another example is the 5 expansion of Bender 
and collaborators jl 1 j . Applied to QED, it can be regarded as an expansion in the number of 
identically charged electron species whose physical value is 6 = 1. Yet an expansion in powers 
of 6 is accurate. It has also been successfully applied to a variety of other non-linear equations. 

Another possible expansion parameter is the inverse number of space-time dimensions 1/d. 
However, we do not yet know of any useful formulation of the d — > oo limit of large A^ Yang-Mills 
theory that is a simplification. This is again motivated by atomic physics, where the d ^ oo 
limit in the zero angular momentum sector is a non-relativistic 0(d) vector model for position 
vectors of electrons. This provides a spectacularly good approximation to the binding energies 
of many-electron atoms in a 1/d expansion, as shown by Herschbach and collaborators jl2j. 

1.2 Loop Equations of Large- Matrix Models 

A primary aim in the study of a Euclidean large- A^ multi- matrix model is to determine its fac- 
torized correlations. They satisfy quantum corrected equations of motion, which are factorized 
Schwinger-Dyson or loop equations (LE). We formulate these in a way that makes manifest some 
algebraic and differential structures they share with the Makeenko-Migdal equations of A^ = oo 
Yang-Mills theoryj5,^^. In particular, they are not differential equations, due to a mismatch 
between the differential and product structures. Though infinite in number and quadratically 
non-linear, we show that they have a hierarchical structure whereby the highest rank correla- 
tions in any equation only appear linearly. However, we show they are underdetermined in many 
interesting cases. We identify additional equations which a naive passage to the large A^ limit 
misses. They are conditions implied by invariance of matrix integrals for correlations, under 
transformations leaving both action and measure invariant, possibly up to 1/A^^ corrections 
(eg. BRST transformations). However, the additional equations are not implemented, so the 
underdeterminacy of the loop equations is not satisfactorily resolved. On the other hand, we 
exploit the algebraic and differential structures to propose an approximation scheme for a class 
of A-(multi)-matrix models motivated by the Lagrangian of Yang-Mills theory. 
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The primary virtue of the scheme is that at zeroth order, it turns the non-hnear loop equations 
into hnear PDEs. Prominent in this class of models are those whose action is a linear sum of 

1 27«- 1 

5g = - tr C'^AiAj, Scs = — C'^^MAj^A^] k Sym = - — tr [^,, A,][Ak, Ai]g'^giK (2) 

In the first two cases, we allow Ai to denote either gluon(hermitian complex) or ghost (grassmann) 
matrices^. Though they arise from terms with 2, 1 and derivatives in the Yang-Mills action, 
these matrix models may be called Gaussian, Chern-Simons and Yang-Mills models since they 
also include the zero momentum limits of the corresponding field theories. The indices i,j,k,l 
are short for position and polarization quantum numbers, while color indices are suppressed. It 
may be possible to fruitfully think of Yang- Mills theory as a grand limiting case of such matrix 
models for appropriate integral kernels C^^,C^^^ and g^^ when the indices become continuous. 
Matrix models and field theories of this type also arise in dimensional reductions of Yang-Mills 
theory to 2 or fewer space-time dimensions. Here we consider bosonic matrix models, the 
extension of our results to models with ghost matrices will be treated in jj^. 

Summary of results and organization: In section l^lD we obtain the large- loop equations^ 
\iJ\S'^^Gji = 6^^^^"^ G i-^G I2 for gluon correlations Gj = {j^ tr Aj) of a hermitian multi-matrix 
model with action tr S{A) = tr Aj . In section 12.21 we show that the loop equations are 
underdetermined in some interesting cases, though they determine infinitely many higher rank 
correlation in terms of lower rank correlations. In section 12.31 we obtain additional equations 
associated with symmetries of both measure and action, which are easily overlooked in passing 
to the large- limit. In section [2.41 the loop equations are reformulated in terms of the series 
G{^) = Gj^^ , where are non-commuting sources: 

J2{n + l)S^^-^-W,,^...D,M0 = G{0eG{0 or S'G{0 = GiOCGiO- (3) 

n>0 

The linear term (variation of action) is written in terms of left annihilation operators Di . The 
quadratic term in gluon correlations involves the concatenation product. It is the variation of the 
matrix model measure and is universal, independent of the action. However, left annihilation 
does not satisfy the Leibnitz rule with respect to concatenation, and to make things worse, 
concatenation is non-commutative. Due to this mismatch, the loop equations are not differential 
equations in the ordinary sense. On the other hand, there is another natural product between 
gluon correlations, the shuffle product (section l2. 5 j) . which arises from the expectation value of 
point-wise products of Wilson loops. It turns out that left annihilation is a derivation of the 
shuffle product. Moreover, there is a democratic version of left annihilation, full annihilation, 
that is a derivation of concatenation (section 12. 6|) . Furthermore, concatenation and shuffle 
combine to form a bialgebra (appendices IbI and IH]) . 

These algebraic and differential structures along with ideas from deformation quantization 
suggest a possible approximation scheme for the loop equations. The idea is to remedy the 
above mismatch by expanding the non-commutative concatenation product in a series around 
the commutative shuffle product so that at zeroth order, concatenation is replaced by shuffle 
and the loop equations become quadratically non-linear inhomogeneous PDEs in an infinite 

^We assume there are an equal number of ghost and anti-ghost matrices in each term, as in Yang-Mills theory. 
^Small letters i denote single indices, capitals denote multi-indices / — i\i2 ■ ■ ■ in and |/| denotes the number 
of indices in a multi-index. Repeated upper and lower indices are summed. S'j is 1 \i I — J and zero otherwise. 
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dimensional space spanned by words in A letters. Thus, the approximation scheme involves the 
introduction of a deformation parameter controlling the amount by which the loop equations 
for gluon and ghost correlations fail to be partial differential equations. The physical value of 
our dimensionless expansion parameter q is 1. 

A further remarkable simplification occurs in models whose action is such that 5* is a 
derivation of the shuffle product. These are models in which 5* is a linear combination of iterated 
commutators of Di and include the zero-momentum Gaussian, Chern-Simons and Yang-Mills 
models as well as their field theoretic counterparts as examples f section l2.7j) . In these cases, the 
passage from G{^) to its shuffle-reciprocal -F(^) = Fj^^ turns the non-linear PDEs into a system 
of linear equations for the Fj (section 14. ip . We obtain an explicit formula for Gj in terms of Fj 
so that once the linear equations are solved, the 0{q^) gluon correlations can be obtained. This 
is illustrated for the zero- momentum Gaussian (section |4.1.1|1 . Chern-Simons fsection l4.1.2|) and 
Yang- Mills (section I4.1.3|) multi-matrix models. For the Gaussian, the linear equations have a 
unique solution which provides a first approximation to the exact large correlations. But for 
the other examples, the equations are underdetermined just as the original loop equations were 
and we exhibit infinite classes of solutions. It remains to find and implemented the additional 
constraints on correlations, such as those associated to symmetries of action and measure. 

In section 14.31 we take the first steps to extend the approximation scheme beyond zeroth 
order. This requires us to find an expansion for concatenation around the shuffle product. Such 
a formula would be loosely analogous to the associative * -product expressions of deformation 
quantization. We obtain two partial results in this direction. First, we find a one parameter 
family of associative g-products that interpolates between commutative shuffle {q = 0) and 
non-commutative concatenation (q = 1). Moreover, by taking q to be infinitesimal, we obtain 
a Poisson bracket on the shuffle algebra. 

In sections 14.21 and 14.3.21 we briefly investigate another approximation scheme for the loop 
equations that involves expanding the left annihilation around full annihilation, holding the 
concatenation product flxed. Though similar in spirit to the main approximation scheme of the 
paper, it has the potential to give a complementary estimate for correlations as shown by its 
application to 1 -matrix models. 

Section EJ is devoted to 1-matrix models. In this case, both concatenation and shuffle are 
commutative, and an explicit 'star product' formula is obtained for the expansion of the former 
around the latter (section l3.2() . In section an expansion for the left annihilation as a series in 
powers of full annihilation is obtained. These lead to two different approximation methods for 
the 1-matrix loop equations, involving either a deformation of the product or the annihilation 
operator. Both schemes are applied to the Gaussian (section 13. 4() . which is the only 1-matrix 
model for which 5* has the derivation property. While deforming the product overestimates 
correlations, deforming the annihilation operator underestimates them. 

Background on Literature: There are several complementary approaches to the loop equa- 
tions of matrix models. First, they are formulated in different ways: resolvents of matrices, 
gluon correlations, planar diagrams, Wilson loops etc. Different approaches to multi-matrix 
models can be broadly categorized by the mathematical structures that play a significant role. 
A major portion of the literature (eg. jl41 1151 ITHl I17j ^ is devoted to exact solutions for certain 
observables of specific (e.g. 1-, 2- and chain-type) matrix models, their multi-cut solutions and 
summing their 1/A^ expansion. This involves connections to integrable systems, algebraic geom- 
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etry and conformal field theory. Another approach exploits the connections to non-commutative 
probability theory (eg. jl81 1191 12UI I21j ) . Yet another point of view seeks to exploit a hidden 
BRST symmetry [221 • A cohomological interpretation of the loop equations and a variational 
principle for them was presented in [20] . The viewpoint in this paper is distinguished by its use 
of algebraic and differential structures and connections to deformation quantization. Its physics 
roots lie in the early work of Makeenko and Migdal|5i iGj, Cvitanovic et. al. [231 01] . loop space 
formalism for gauge theories |25 | I26 | [27] 128] . and the more recent investigations of Rajeev and 
coworkers [2^1 ISOllS HOI ISllS] • Some structures used in our constructions (eg. shuffle products 
and their deformations) appear in the mathematics literature on calculus of loop space due to 
Chen |32j . the theory of free Lie algebras ^3 and the deformation theory of (Hopf) algebras 
|34| 135] . A feature of the present work is that we do not make any a priori restriction to a 
subclass of correlations (eg. 'mixed' or 'unmixed') as is often assumed in the literature. 



2 Algebraic structure of loop equations of multi-matrix models 
2.1 Factorized loop equations for gluon correlation tensors 

We begin by obtaining the loop equations of a bosonic multi-matrix model in terms of gluon 
correlation tensors. This is convenient to study their algebraic structures and permits treatment 
of all factorized N = oo correlations without restriction. Consider a Euclidean A-matrix model 
with polynomial action tr S{A) = tr S"^ A j . Let = tr Aj denote the 'loop' variable. The 
partition function and gluon correlations are 

Z = J n.dvl.e-^ ^(^) and {<^k, ■ ■ ■ '^kJ = | / n^d^.e"^ '^^^<^k, ■■■<^k„. (4) 

Gk = ^^^n—>oo{^k) are the gluon correlations of interest in the large- limit. Here Ai = 
A\, 1 < i < a are N x N hermitian matrices. The tensors S are the 'coupling tensors' defining 
the theory. Due to the trace, the only part of 5^ that contributes is its cyclic projection, so 
assume that are cyclically symmetric, 5"^* = S^^ for all Gluon correlation tensors Gj 
are also cyclically symmetric. Additionally, assume are chosen such that (S^)* = where 
7 is the word with indices reversed^. This, along with hermiticity of Ai ensures that tr S{A) 
is real. In turn, this implies that G*j = Gj. To see this, recall that for any complex matrix M, 
( tr M)* = tr M"^ and apply this to M = Aj and use hermiticity of Ai . For the Gaussian, all 
= except S^^ which may be taken as a (positive) real-symmetric matrix. 

The Schwinger-Dyson equations(SDE) are constraints on {^Ki '''^Kn) implied by invari- 
ance of the matrix integral under an infinitesimal (but non-linear) change of integration variable 

[Ai]'^ [A'i]'^ = [Ai]'^ + vf[Ai]'^, where Vi are infinitesimal real parameters. (5) 

Under this change of variable, the infinitesimal changes in , the action and the measure are 

det(llll) = 1 + aV5{-^^c^,,<I>,,. (6) 



^This is satisfied by examples such as the Gaussian, Yang-Mills and Chern-Simons theories, see Sec 12.71 
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Invariance of {^Ki ■ ■ ■ ^Kn) to linear order in v( implies the SDE^ 

/ n 

viS'^^'^H'^j.u,) = vi6'j^''^'^i,<^j,) + ^Y.^if'"i'^L,iM,), V K, and n = 0,l,2,... (7) 

p=i 

So far we have not made any approximation. In the large limit, expectation values of U{N) 
invariants factorize ($/^<l>/2) = {^Ii){^i2) J.- Naively, the leading factorized Schwinger-Dyson 
or loop equations (LE), which are a closed system for G/, are 

v',S'^'-'^Gj,ij,=vi5'/'''Gi,Gi, yv (8) 

These infinitesimal changes of variable are associated to vector fields = v(L^j whose action 
on Gj is given by L\Gj = Sj^^'^^G j^ij^ ■ In particular, choosing the components of the vector 
fields v( to be non- vanishing only for a single (i, /) , we get the loop equations 

S^''-^'Gj,ij,=6'/'''Gi,Gi, V/, i. (9) 

Using cyclicity of and G/ we get 

\iJ\ S-^'Gji = 6j^'^^Gi,Gi,, V /, i. (10) 

LE 1)10(1 relate a changes in (expectation values of) action and measure under the action of L\. 
However, there may be vector fields (i.e. choices of v( ) for which both sides of |H1 vanish^ . In 
that case, the leading equation in the large limit is different from (see section IT^ . 

We seek solutions to El among cyclic symmetric tensors Gj satisfying G} = Gj and G0 = 
Go = 1 , where is the empty string. Note that the LE may make sense even when the matrix 
integrals don't seem to converge, as for a cubic action. When analogues of (|1U() are formulated 
for Wilson loops in a gauge theory [S], they are called Makeenko-Migdal equations (notice the 
resemblance between (fTUI) and dJ) 

'^^^^^(^) = ^ / dyj^^^x - y)W{Gy,)W{G^y). (11) 

2.2 Underdetermined nature of loop equations and examples 

Given an action S{A) , G/ are uniquely defined by (HJ provided the integrals converge. As 
examples below show, the large- LE (|10() do not determine G/ uniquely in general. In section 
1231 we obtain additional large- SDE involving Gj that were not accounted for in the passage 
from ^ to (|10() . But even these may not be sufficient to fix the G/. 

Consider first A = 1 matrix models whose LE are got by restricting (|10|) to a single matrix. 
Suppose tr S{A) = tr X^I^i ^lA^ is an m^^ order polynomial, then if G^ = {^A'') 

m 

Y,lSiGk+i= J2 GrG,, for k = -1,0,1,- ■■■ (12) 

1=1 ''i'5>0, r+s=k 

^Sometimes called Virasoro constraints in string models or Ward identities. Ward identities seems more 
appropriate to the special case where the change of integration variable was a gauge or BRST transformation. 
^Note that this may happen even if there is no (z, J) for which both sides of llOl vanish. 
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The LE listed sequentially are 



k = -l: Si + 2S2Gi + --- + mSmGm-i = ^, 

k = Q: SiGi + 2S2G2 + ■■■ + mSmGm = 1, 

k = l: SiG2 + 2S2G^ + --- + mSmGm+i = 2Gi, 

k = 2: SiG3 + 2S2G4 + --- + mSmG^+2 = 2G2 + Gl ... (13) 

We see that in the k^^ equation, the highest rank correlation G^j^}« appears linearly [Sm 7^ 0) 
and may be determined in terms of lower rank correlations. For a Gaussian (m = 2) 1)12(1 
determine all moments. More generally, the LE determine higher moments Gm-i, Gm, Gm+i, ■ ■ ■ 
in terms of m — 2 undetermined lower moments Gi, . . . Gm-2 ■ However, among Gi, • • • Gm-2 , 
the odd ones must vanish if the action is even. Observe that this is associated with the [A]"^ 1— > 
— [A]^ symmetry of an even action and of the measure if — > 00 through even values. Such 
transformations provide additional equations missed out by the LE. 

For multi-matrix models, suppose S{A) is an m*^ order polynomial, i.e S'^ = if \ J\ > m 
and 3 J with \ J\ = m such that S"-^ / 0. Then the loop equation \iJ\S-^^G ji = Sj^'^'Gi^^Gi^ for 
any fixed / and i involves correlations with highest rank (|/| + m — 1) only linearly. Of course, 
there are several correlations with a given rank and several equations for fixed |/|. If all Gk 
up to \K\ < r are known, we have a system of inhomogeneous linear equations for correlations 
of rank r + 1. For the Gaussian tr S{A) = ^ tr G^^AiAj, these are just recursion relations 

Gij = GijS^^-^^^Gi^Gi^ where GijC^'' = 5^ . Their unique solution for all correlations is given by 
the planar version of Wick's theorem, which is a sum over all non-crossing partitions of il into 
pairs. But for many interesting cubic and higher order actions, the LE are underdetermined 
even by comparison with 1-matrix models. Not only are Gk for \K\ < m — 2 left undetermined, 
many higher rank correlations are also not determined in terms of them. Consider two examples: 
a quartic 2 -matrix model and the Chern-Simons 3 -matrix model. 



2.2.1 Quartic 2-Matrix Model 

Suppose tr S{A) = tr [CA1A2 + ^{Af + ^2)]- The matrix integrals converge and the cyclic 
couphng tensors are S^^^^ = 52222 _ g ^^^^ _ ^21 _ c ^ r^j^g £qj, ^^^^ j 

cG2i + gGuii = S'/^''Gi,Gi, and cGu + 9G2221 = S'/^''Gi,Gi,. (14) 

Since the action is an m = 4*'^ order polynomial, the LE do not fix Gi, Gij . They determine an 
infinite number of higher rank correlations in terms of these, but also leave an infinite number 
undetermined. For / = the two LE give Gm = ~|G2 and G222 = The other rank-3 

correlations Gii2,Gi22 are left undetermined. For I = ii, the LE determine 4 of 6 correlations 
leaving G1122 and G1212 undetermined: 

1 c c 

Giiii = G2222 = -(1 — CG12), G1112 = — G22, G1222 = — Gil. (15) 
9 9 9 

For I = iii2, the LE are 

cG2i^i2 +5^1111112 = ^h'^h ~^ ^iiGi2 and cGij^jj +5^222*112 = ^iiGh + ^f^Gi^. (16) 
They determine 6 of the 8 rank- 5 correlations in terms of lower rank ones 
1 

Giiiii = -(2Gi — CG112), G11112 = -(G2 — CG122), G11122 = ^Gi, 
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^22222 — -(26*2 — cGi22) ^12222 — 0(^1 — CG112), ^11222 — -5-^2, (17) 

while leaving G12121 and G21212 undetermined. In this manner, by choosing longer words /, we 
can fix an infinite number of higher rank correlations in terms of lower rank ones, but at each 
step a few correlations remain undetermined. The number of undetermined correlators may be 
significantly reduced by the Ai ^ A2 symmetry of S{A) which implies Gj = Gj if / can be 
obtained from J by 1 <-> 2 and a cyclic permutation. Notice that this is also a symmetry of the 
integration measure. The same applies to the change of variables Ai i— > —Ai,A2 1— > —A2. 

2.2.2 Chern-Simons Model 

The LE of the CS model tr S{A) = ^e*-''^ tr AiAjAk are 

2iKe'^''Gijk = S'/'''Gi,Gi,. (18) 

They leave rank-1 correlations Gi undetermined (m = 3). For |/| = and arbitrary i, the 
LE are e^^^Gjk = which do not give any constraints not already implied by cyclic symmetry 
of Gjk- Thus G12, Gi3, G23, Gil, G22, are all left undetermined. For |/| = 1 with arbitrary 
I = ii and i, the LE are 2iKe^^^Gjki^ = 6l_^ . From 9 possible (complex) equations we get only 
1 independent condition after accounting for cyclicity and hermiticity: the imaginary part of 

G123 — Gi32 = (19) 

ZlK 

This allows us to fix only one parameter in the c(3, A = 3) = 11 dimensional space of 3^^^^ rank 
cyclic hermitian tensors (see appendix^. For I = iii2 and i arbitrary, the LE are 

2iKe'^''Gi,i,,k = 4 Gi, + 4 Gi, . (20) 

Of the 27 possible equations, there are actually only 9 independent ones that do not follow 
from cyclicity^. Three 'homogeneous' ones G1212 = G1122, G1313 = G1133 , G2323 = G2233 and 
six 'inhomogeneous' ones 

2iK(Gii23 - G1213) = Gi, 2iK(Gi2i3 - G1132) = Gi 

2iK(Gi223 — G1232) = G2, 2iK(Gi232 — G1322) = G2 

2iK(Gi323 - G1332) = G3, 2m(Gi233 - G1323) = G3. (21) 

Nevertheless, these conditions are not enough to fix the c(4, A = 3) = 24 independent cyclic and 
hermitian 4**^ -rank tensors (see appendix^. This underdetermined nature of the LE persists 
for correlations of higher rank. Notice also that by Ai ^ A2 ^ A^ ^ Ai symmetry of the 
action and measure, we have Gi = G2 = G3 etc, but this is not a consequence of the LE and 
still leaves the common value of these undetermined. 

®The fact that many of the loop equations are not independent of each other indicates there are vector fields 
for which both sides of|H|vanish identically. 
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2.3 Additional equations for gluon correlations 



Are there more equations satisfied by G/ that wih lessen the underdeterminacy of the LE? 
In going from finite- SDE ^ to large- LE H1U|1 . we overlooked the possibility that both 
LHS and RHS of Q may vanish for some v . In other words, Ai ^ Ai -\- vl Aj may leave the 
(factorized expectation value of) action and measure simultaneously invariant at leading order 
as N ^ oo. For such v( the 0{N^) terms in identically vanish and the 0(1/A^) terms 
constitute the leading large- A^ SDE. Denote 

G?^ Gf^ Gf\ Gf, 

Then the ©(l/A^) terms in ((TJ b ecome 

n 

vis'^^'^Gfjj^ = vi5\'^'^Gflj^ + vi bX'^'GL.iM, V V, K, and n = 1, 2, . . . (23) 

p=i 

Unfortunately, (|23() involve not just the Gi but also 1/A^ corrections to single and double- 
trace correlations. Thus, an attempt to ameliorate the underdetermined nature of the LE seems 
to open a new can of worms. However, in keeping with the spirit of the large- A^ limit as an 
approximation where we retain only the leading large- A^ contribution to all quantities, it seems 

(2) 

reasonable to ignore the G... terms and consider 

n n n 

J2 vis'i.f'^GL.iM, = ^^Y. ^iL\GK, = ^^Y. ^'"Gk, = (24) 
p=i p=i p=i 

At first, these equations seem universal, they do not involve the coupling tensors at all! 
However, for generic v, these are 1/A^^ contributions to the SDE and should be ignored in the 
large- A limit. But if are such that both RHS and LHS of Q vanish identically, then these 
become the leading large- A SDE. Thus, these equations are not universal, since they must be 
enforced only for those vf for which the leading change in action and measure vanish identically. 
To summarize, the additional equations are 

n 

^L^Gxp = V Ki,---,Kn and n = 1,2,3... 

p=i 

and all vj such that v{ S-^'^'^^Gj.ij, = vj5j''^^Gi,Gi^ = 0. (25) 

Are there any such additional equations? This is related to whether there are any transforma- 
tions that leave both action and measure invariant at leading order as A ^ oo. We exhibited 
several such discrete transformations in sections 12.21 12.2.11 and 12.2.21 BRST transformations of 
gauge fixed Yang-Mills theory are also of this sort and lead to Ward or Slavnov- Taylor iden- 
tities. Are the LE (|in|) consistent with the additional equations (|25|) ? This would vindicate 
our throwing away the subleading G... terms in (|23|1 . If so, do the LE (|l()jl together with (|25|) 
determine the Gj , or do we need yet more conditions? We postpone investigation of these very 
interesting issues and focus on the LE in the rest of this paper. 
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2.4 Loop equation in terms of left annihilation and concatenation 

Define the generating series of gluon correlations by the formal sum G{^) = Gj^^ . Here, 1 < 
i < A are non-commuting variables that can be thought of as sources, and = . . . . 

If they did commute, the generating series would only contain information about the symmetric 
correlations. But since are not symmetric in general (only cyclically symmetric), there 

is no relation between ^'^-^ and ^^^"^ . Define the concatenation product cone by 

^i^J = or = F/GjC" ^ iFG)K = S'j/FiGj. (26) 

For example'', 

{FG)o = FoGo; (FG), = F,Go + FqG,; {FG),, = F^G,, + F,G, + F,,Gq- etc. (27) 

In terms of cone, the RHS of ^ becomes 6j^'^^Gi,Gi^ = [G{OCtG{0]i ■ Also define left 
annihilation^ 

Dj.f = 5*if2-»". (28) 
Dj eliminates the left most source if ii = j and returns zero otherwise. In terms of coefficients, 

[D,G]j = G,j, [D,„ . . . D,,G]j = Gj,...,„j, (29) 
so that Gjj = [DjG]i . The LE (|10|) . one for each i, can be written as 

Y^in + l)S^^-^"^D,„...DjMO = GiOeGiO or 5^G(e) = G(Of G(e). (30) 

n>0 

We used cyclicity of S^,Gj in deriving this. Thus, the LE involve left annihilation and cone 
product. The LHS of (PU)) defines the action dependent operator 

S^ = J2in + l)S''-'"'Dj„...D,,. (31) 

n>0 

At first glance, the LE (|3()|) look like quadratically non-linear PDEs whose order is one less than 
that of the action polynomial. However, concatenation in the universal term on the RHS is 
non-commutative since sources ^* do not commute. Further, left annihilation does not satisfy 
the Leibnitz rule with respect to concatenation, i.e. Dj are not derivations of cone. This 
'mismatch-match' between product and annihilation make the LE difficult to solve. It turns out 
there is another natural product between gluon correlation tensors, the shuffle product, with 
respect to which left annihilation satisfies the Leibnitz rule. We try to exploit the interplay 
between cone, shuffle and their derivations to find an approximation method to solve the LE. 

2.5 Shuffle multiplication from products of Wilson loop expectation values 

Here we obtain the shuffle product of gluon correlations induced by expectation values of prod- 
ucts of Wilson loops. The expectation value of the Wilson loop F{'y) is a complex- valued 

''Note that concatenation of cyclically symmetric tensors is not cyclically symmetric in general. 
*Left annihilation does not preserve cyclic symmetry of tensors in general. 
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gauge-invariant function on the space of loops 'y : M , where M is space-time. If Ai^{x) 

denotes the components of a gauge field 1-form valued in the Lie algebra of hermitian matrices, 
we define the path ordered exponent 



^(7) = -^trT'exp 



Jo 



i I A,Ax)^^ds 
as 



(32) 



Parameterized loops on M are denoted x''{s). Wilson loops are typical functions on loop-space 
and their expectation values can be expanded in iterated integrals of gluon correlations 



iHl)) = E^™/ {-tTA,,{x{s,))...A,MSm)))-T-----j — ds,...ds^ 

^0 •'0<si< - <sm.<l dSi dSm 

= y^i"' Fvi-ym{x{si),--- ,x{sm))—, 3 dsi---dsm (33) 

^0 ^0<si<-<s^<l dsi dSm 



where the gluon correlation tensors associated to -^(7) are 

F^,...^Msi),---,x{sm)) = {^tiA^^{x{si))---A^Msm)))- (34) 

The point- wise commutative product of functions on loop-space is defined as (FG) (7) = F{'y)G{'j) . 
Taking expectation- values and working in the large- iV limit, where correlations factorize, we get 

((FG)(7)) = (F(7)G(7)) = (i^(7))(G(7)) + 0{^). (35) 

We may expand the LHS in correlation functions associated to the Wilson loop (FG) (7) . We 
call these {F o G)pj^...pp(x{ui) ■ ■ ■ x{up)) . They are defined as 

((FG)(7)) = E / {Fo G),,..,^(x(ui) • ■ • x{u,))'^ ■ ■ ■ ^dm ■ ■ ■ du,. (36) 

^ ./o<ui<-<up<i dui dUp 

Meanwhile, the expansion of the RHS reads 

CO „ 

(F(7))(G(7)) = E ^"^^ /o<.i<....™<ii^:.i.....(^(si),---,a;(s^))G^i...,.„(^(ii),---,a;(in)) 



m,n=0 



dx^^ dx""" dxf^^ dx^" , , , , 
X — Ti—dsi ■ ■ ■ dsmdti ■ ■ ■ dtn- (37) 

dsi dSm dti dtn 

To make this look like the expansion of the LHS, we collect terms with a common sum n+m = p 
and then sum from p = to 00. Moreover, we must relabel the i^'s and //'s as p's and the s's 
and t^s as u 's. We must allow every possible relabeling that preserves the order among the s 's 
and t's. When this is done, we read off the relation between the gluon correlations associated 
to the Wilson loop {FG){j) and those associated to -^(7) and 6^(7) 

{FoG)p,...p^{x{ui)---x{up)) = E E ^P,-i(i)-P,-i(^)(a^K-i(i))r--,a;(n^-i(^))) 

m+n=pa an (m,n) 

shuffle 

^'^Pa-Hm+iy-P.-Hm+nM^'r-Hm+l)), ' ' ^K"! (m+n) )) (38) 

An (m, n) shuffle is a permutation of m -|- n letters {1,2, - ■ ■ ,m + n) such that 

cr~^(l) < • • • < o-~^(m) and a~^{m + !)<■■■< a~'^{m + n). (39) 
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For brevity, we combine the Lorentz and space-time indices into a single index i, then 

^ ' -L / J / J (J ^(1)' cr (m) (J -"-(m + l)' a ^ (m + n) ^ ^ 

m+n=p a an {m,n) shuffle 

The RHS is called the shuffle product (s/i). It is commutative. A compact notation for sh is 

{FoG)i= ^jGk. (41) 

I=JUK 

The condition I = J UK means that J and K are complementary order-preserving sub-words 
of /. The operation J U K is a riffle-shuffle of two card packs J and K . Some examples are 

[F o G]i = FiGo + FoGi] [F o G]ij = FijGo + FiGj + FjGi + FoGij; 

[F o G]ijk = FijkGo + FijGk + Fii^Gj + FjkGi 

+FiGjk + FjGik + FkGij + FoGijk] 

[F o G]ijki = FijkiGo + FijkGi + FijiGk + Fi^iGj + Fj^iGi 

+FijGki + FikGji + FiiGjk + FjkGii + FjiGik + F^iGij 

+FiGjki + FjGiki + FkGiji + FiGijk + FoGijki- (42) 

We notice two properties of sh. If Fj and Gj are cyclically symmetric for all / and J, then so 
is {FoG)k for all K . To see why this is true in general, observe that {FoG)k is the expectation 
value of the trace of a product of gluon fields, and the trace makes it cyclically symmetric. Thus 
sh preserves cyclicity of tensors. Moreover, we notice that if Fj and Gj satisfy the hermiticity 
properties Fj = Fj, Gj = Gj for all /, J, then so does their shuffle product 

(F o G)} = {Fo G)j V /. (43) 

This is a reflection of the relations^ ^(7)* = ^(7) and {FG)*{-f) = F*{'y)G*{-f) = {FG){j) 
when the path-ordered exponential is expanded out in iterated integrals. 

The shuffle product allows us to reduce manipulations in the commutative algebra of func- 
tions on the infinite dimensional space Loop{M) to operations on tensors on the finite dimen- 
sional space M . More precisely, start with a manifold M, and denote the space of 1 -forms on 
M by A^(M). Then consider the tensor algebra T on A^(M). The shuffle algebra is 

Sh{M) =T{A^{M)). (44) 

The shuffle algebra is a replacement for the algebra of functions on Loop{M) . Let , , . . . 
be a basis for A^(M) (think of these as (ix*^, . . .), then an element of the shuffle algebra is 

G = J2 Gn-iX' ® • • • = E G',,..,„r^-"", (45) 

n n 

and is to be regarded as a function on Loop{M) . A specific collection of gluon correlations 
{Gii...i^}^Q can encode the information contained in the expectation value of a specific function 
G(7) on Loop{M)^^ . The shuffle product of basis elements is 

f o = f J' + ^J'*; f J' o = f J''^ + f '^J' + ^''^^ (46) 
^7 is the loop 7 with opposite orientation. 

^"The map is not 1-1 since gluon correlations are not gauge invariant in general, unlike Wilson loops. A way 
to deal with this is to introduce ghosts. When LE are formulated in terms of Wilson loops, gauge fixing and 
ghost contributions cancel out jTj. But this is not the case if we work with correlation tensors. Extension of this 
formalism to include ghosts in matrix models will be treated in ^33- 
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and in general 

^h-ip Q ^ip+i-ip+q ^ ^ ^V(i)-V(p+g) Qj. o = . (47) 

a a (p,q) shuffle 

To summarize, we have shown that the commutative point-wise product of Wilson loops induces 
the commutative, cyclicity and hermiticity preserving shuffle product of gluon correlations^^ . 



2.6 Derivations of shuffle and concatenation products 

Concatenation and shuffle combine to define a pair of dual bialgebras on the vector space 
span(^^) (see appendices IbI and IH]). Derivations of concatenation and shuffle play a central 
role in this paper. Recall that the LE H30|) involved left annihilation Di defined in (|28() . We 
show that Di is a derivation of sh i.e. it satisfies the Leibnitz rule 

Di{FoG) = {D,F)oG + Fo{DiG). (48) 

The proof is by explicit calculation [Di{F o G)\i = [F o G\ii = J2iiui2=ii ^h^h ■ Now either 
i £ Ii or i £ I2, so 

mFoG)]i = FihGi,+ E FhGu,= E [D^F]I,GI,+ ^/ilAG]/, 

/lU/2=/ hUl2=I IlUl2=I hUl2=I 

= [{DiF)oG]i + [Fo{DiG)]i. (49) 
Full annihilation^^ Dj is a democratic version of left annihilation. It is defined as 

^3^' = ^hjh^'''' and [T>jF]i = 5\^'^Fi,ji,. (50) 
Dj does not preserve cyclic symmetry of tensors. However, Dj is a derivation of cone, 

T>j{FG) = (DjF)G + F{T>jG). (51) 
To see this, begin with the LHS [Uj{FG)]i = 6j'^\FG)i,ji^ , 

[D,(FG)], = S'/'^5f^]l^FK,GK2 = <5^^^^^FiuX2GL3 + 5f i^^^^F^.G^.^Xa- (52) 
On the other hand, 

[(D,F)G], = 6'j^''{B,F)j,Gj, = 5]^''6il'^Fj,,j,Gi, = S^^'^^'^^Fl^.l^Gl,. (53) 

Thus 

[(D,F)G], + [F(D,G)], = 5^^^^^Fi,,x,Gz.3 + ^^^^^Fz.,Gi,,x3 = [D,(FG)],. (54) 

The commutator of derivations is a derivation irrespective of whether the product is commutative 
or not. This is analogous to the Lie bracket of vector fields being a vector field on a manifold. 
For example, merely using the fact that Di is a derivation of s/i = o, it is easy to show that 

[D,,D,]{FoG) = {[D,,D,]F)oG + Fo{[Di,Dj]G). (55) 

^^This construction generalizes to differential forms on Loop{M) , but we do not use it in this paper. 

^'^ Not the cyclic gradient. The cyclic gradient is 5i£,' — Si-^n^^.^^^^ and is not a derivation of concatenation. 
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It follows that iterated commutators of derivations (e.g. [Di, [Dj,Dk]]) are also derivations. On 
the other hand, products of left annihilation operators are not derivations of the shuffle algebra. 
For e.g. DiDj = Dij is not a derivation of sh. This is analogous to the product of vector 
fields not being a vector field. Furthermore, left annihilation operators with a single index Di 
do not form a Lie algebra by themselves. The commutator [Di,Dj] = Dij — Dji is not a linear 
combination of Df. 's. However, by construction, the vector space spanned by the set of all 
iterated commutators of left annihilation operators Di,[Di,Dj],[Di,[Dj,Dj^]], - ■ ■ forms a Lie 
algebra, the Lie algebra of derivations of the shuffle product. This is the free Lie algebra. It is 
analogous to the Lie algebra of left invariant vector fields on a Lie group. Here, the role of the 
Lie group is played by the free group on A generators. 

2.7 Derivation property of terms in Yang-Mills action 

The action-dependent linear term 5*G'(^) in the LE (|3U() is a sum of products of left annihilation 
operators 5* = J2n>o{''^ + ^)^''^ ''"^^jn ' ' ' ^ji ■ Suppose coupling tensors are such that 5* 
is a linear sum of iterated commutators of left annihilation operators, 

= C'W, + [D, , Dk] + &^'''[[Dj ,Dk],Di] + --- (56) 

Then 5* is a derivation of shuffle. Of what practical use is this property? The LE (|3U() are 
quadratically non-linear in cone, but involve left annihilation, which is a derivation of sh. In 
section^ we introduce an approximation scheme where cone is expanded around sh. The main 
simplification for matrix models having the derivation property is that their LE can be turned 
into (an infinite system of) linear PDEs at 0**^ order in this approximation. This is not the case 
for matrix models without the derivation property. 

Among 1 -matrix models, the only one with this property is the Gaussian tr S{A) = 2^ tr ^4^ 
for which S = ^D . For A = 1, there is only one left annihilation operator, and all its iterated 
commutators vanish. Multi-matrix models provide non-trivial examples. It is remarkable that 
the gluonic terms in the Yang- Mills action quadratic in momentum, linear in momentum and 
independent of momentum each separately has this derivation property^^. These terms can be 
written as tr C'^AiAj, tr C^^^Ai[Aj,Ak] and tr [Ai, Aj][Ak, Aijg^'' g^'- for appropriate tensors 
C^i ^ C^i^ ^ g^i . Moreover, the zero momentum limits of the Gaussian, Chern-Simons and Yang- 
Mills matrix field theories all have this derivation property. They correspond to the simplest 
non-vanishing choices for the tensors C*-' , C"^^^ , (jijkl (j5(jj) . In fact, this property also extends 
to the corresponding matrix field theories but we do not address that here. 

Gaussian: The Gaussian multi-matrix model tr S{A) = ^ tr C^^ AiAj has real-symmetric 
covariance C*-' = C-'*. S^^ = ^C*-' is cyclically symmetric and also satisfies (S'*-')* = S'-'* so that 
all correlations satisfy G} = Gj. We get 5* = 2S^^Dj = C^^Dj, which is a linear combination 
of left annihilation operators and therefore a derivation of sh. The LE are 

C^W,G{C) = GiOCGiC). (57) 

Chern-Simons: For at least three matrices (A > 3) , the CS type of matrix model has action 
^ tr G^^^Ai[Aj,Ak] where C^^^ is any tensor which is anti-symmetric under interchange of 

^^See for the corresponding property after inclusion of ghosts. 
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any pair of indices. The part of C^^^ that is symmetric under interchange of a pair of indices 
does not contribute on account of antisymmetry of the commutator. The action can also be 
written as tr S{A) = ^ tr C'^'^AiAjAk where C'^^ = C'^^ - C^^^ . The particular case of zero 
momentum 3d CS gauge theory results from the choice A = 3, C^^^ = e^^^ (the Levi-Civita 
symbol), and integer- valued coupling constant ^ttk. More importantly, terms in the Yang- 
Mills action linear in momentum are of this form. Irrespective of its field theoretic origin, 
S'^^ = {2iK/3)C'^^ is cychcally symmetric since C'''^ = {-ifC'^''. Moreover, {S'^'')* = S^^^ so 
that G*j = Gj . Now 5* is a linear combination of commutators of left annihilation operators: 

S' = 2iKC'^^DkDj = iK{C'^''DkDj - &^WkDj} = in&^^Dk, Dj] (58) 

and therefore a derivation of sh. The 'Chern-Simons' loop equations are 

iK&^\Dk. Dj]G{i) = G{i)eG{i). (59) 

Yang-Mills: For A > 2, the zero momentum limit of Yang-Mills theory has action (a = g^) 

tr S{A) = tr [Ai,A,][Ak,Ai]g'^gi\ (60) 

where g^^ = g^^ is the inverse metric, it is a real symmetric matrix. The action is rewritten as 

tr SiA) = ^ tr ig^'g^' - g''g^^)A,,u = tr m^g^' - g^'g^^ - <7^^/V.,h] (61) 

so that S''^^'- = -j^{2g^^g^'--g'''-g^'^-g'^^g^^) is cyclically symmetric. Moreover, S'*-''^' = (5"«i«)* = 
gikji follows since g^^ is real symmetric. Then the differential operator 5* = {3 + 1)3'^^^^ DiDj^Dj 

= -^g'^'g^'iDiDkDj - D^DiDj + DiDkD, - DiD^Dj,) = -}-g'^g^\D.j, [D^, A]] (62) 

is a linear combination of iterated commutators of derivations and hence a derivation of the 
shuffle product. The Yang-Mills LE are thus 

- ^sVfl^j, [Dk,Di]]G{i) = G(OrG(e). (63) 

On the other hand, most matrix models do not have this derivation property. For example, 
consider the popular two matrix model tr 5(j4i, A2) = tr [A\ + A2 + 2AiA2\. Here, 
5^ = 2D2 + 4L>f and 5^ = 2Di + 4L>| are not linear combinations of iterated commutators of 
Di and do not define derivations of the shuffle algebra. 



3 Approximation method for one-matrix models 

The LE of a 1 -matrix model ()12|) with m*'^ order polynomial action 

m 

Y,lSiD'-^G{0 = GmG{0- (64) 
1=1 

can be written in terms of left annihilation^^ D . Concatenation, which appears on the RHS is 
the usual product of calculus. But D satisfies the Leibnitz rule with respect to sh, not cone. So 
this is not a differential equation. We develop approximation methods to solve these LE either 
by expanding cone around sh or by expanding D around full annihilation D H50|). which is a 
derivation of cone. Both these turn the LE into linear ODEs at each order of the expansion. 
^''The 1-matrix left annihilation operator, D^" = ^n-i ^-^^ same as the usual derivative of calculus. 
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3.1 Shuffle, concatenation and their derivations 

We give the 1 -matrix versions of cone, sh and their derivations by speciahzation from sections 
12.41 and 12.51 Then we define g -deformed products and derivations that we use to solve the LE 
approximately. Suppose F{^) = J2n>o^n^^ etc. Cone = *i is the usual product of calculus^^, 

n 

^P*,^i = ^P+<i or (F*iG)„ = 5^F,G„_, (65) 

r=0 

while shuffle = *o (previously denoted o) is, 

e *0 ^'=(^1 ^'^'^ or {F *o G)n = E (^^ FrGn-r (66) 

For example ^ *o = . Both are commutative. The notation anticipates *q that interpolates 
between sh {q = 0) and eone {q = 1). We also define 1-matrix analogs of left and full 
annihilation and name them in anticipation of g-annihilation Dg. Left annihilation Dq^^ = 
is the 1-matrix version of Di defined in (|28l) . Dq is a derivation of shuffle 

{DoF)n = Fn+U {Do{F *o G))n = {{DoF) *o G)n + {F *o iDoG))n. (67) 

Full annihilation Di^^ = n^"~^ is the same as the usual derivative of calculus. It is the 1-matrix 
version of Dj defined in (|5fl|) . Di is a derivation of eone, 

[DiF]n = {n + l)Fn+i, {Di{F*iG))n = {{DiF)*iG)n + {F*i{DiG))n. (68) 
This follows from the easily verified formula 

71+1 n n 

in + l)J2 FrGn+l-r = + l)^r+lG„-r + ^(n - r + l)FrGn-r+l. (69) 

r=0 r=0 r=0 



3.2 g-Deformed product 

The g-product interpolates between eone {q = 1) and sh {q = 0)^^ 

(F G)n = J2 \] FrGn-r. (70) 
r=0 V ' 1-q 

It is associative and commutative for < g < 1. The (7 -binomial coefficients or Gauss binomials 
(")^ are polynomials in q with non-negative coefficients. They reduce to unity for q = and to 
the usual binomial coefficients when q = 1. To obtain their properties let yx = qxy . Then 

(x + y)" = X:H ^""V. (71) 

r=0 \^/ q 

^®We also denote ccmc = *i by juxtaposition. 

^^The quantity 1 — g often occurs in formulae, so we call it p = 1 ~ q. 
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The first three Gauss binomials are 



n\ ^ n 



^n^ _ J(l + g2 + g4_^ ^gn-2)(x + g + g""^), if n is even 

(1 + g2 + ^4 _^ . . . + g«-3)(l + g + g2 ^ . . . + ■£ ^ 



The g-Pascal relation is got by multiplying (x + y)"^^ by (x + y) either from the right or left: 

Substituting the first in the second gives 

(:)^ = I^(";')^ fo. 0<.<. (74) 

Iterating, we get 



(75) 



This can also be written as 



n \ \n\ 



The symmetry (")^ = (^"j,) is now manifest, which guarantees commutativity of the g-product 
((701) . Some examples of the g -product are 

(F*,G)o = FoGo; (F *g G)i = FiGq + Fod; (F *g G)2 = FqGs + (1 +p)FiGi + F2G0; 
(F*,G)3 = FoG3 + (l+p + p')(FiG2 + F2Gi) + F3Go; 
(F*,G)4 = FoG4 + (l+p + p2+p3)(FiG3 + F3Gi) 

+{l+p + 2p^ +p^)F2G2 + F^Go. (77) 

We expand the g -binomials around the ordinary binomial coefficients (g = 1) in a Taylor series 

n\ I n\ { r(n — r) oxi x 

> - ' ^ - \ ' p + 0{p^)\. (78) 



\ / (J \ / 1 
Thus *q may be expanded around shuffle *o 

(F*,G)n = {F*oG)n-^J2\]'-Pr{n-r)Gn-r + 

^ r=0 y ) 1 

n 



(F*oG')„-f^ [e*0^0F(e)]r[e*0^0G'(O]n-r + 



9 

^ r=0 V / 1 

(F*gG)(e) = (F*oG)(e)-| e*o(^oF)(e)*oe*o(I)oG)(0 + ---. (79) 
Taking g = 1 , and using commutativity of *o , we get an expansion for cone in terms of sh 

(F *i G) (X) = (F *o G) (e) - 1 e *o C *o (^oF) (C) *o (^oG) (0 + • • • • (80) 
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3.3 g-Deformed annihilation operator 



Recall from section [3.11 that left annihilation [DQF]n = i^^+i and full annihilation [DiF]n 
Fn+i . More generally, let 



{DqF)n =[n+ l]qFn+l 



-qn+1 _ y 






Fn+1 — 


[ q-l \ 





1+q + q^ + ■■■ + q'" 



F 



n+l- 



(81) 



Dq reduces to left and full annihilation for q = and q = 1. However, Dq is not a derivation 
of *q for < q < 1. Fortunately, we don't seem to need that. More importantly, we expand 
Dq around Di in powers of p = 1 — q. Denoting cone reciprocal by usual division of calculus. 



k=l 



k-l 



1 



D^Fii) 



(82) 



3.4 Gaussian one matrix model 



Now we apply this formalism to the simplest of matrix models, the Gaussian 1-matrix model. 
We pick it as it is the only 1-matrix model with the derivation property. We show how expanding 
cone around sh and expanding Dq around Di , are used along with the derivation property to 
turn the non-linear LE into linear ODEs at each order in our approximation schemes. The 
resulting gluon correlations are compared with the exact solution. 

From H12|) . the LE for the Gaussian 1 matrix model with action = ^ tr are 
Z)o^(0 = *i e *i ^(0 or Gn+i = a ^ n = 0,1,2,- •• (83) 

r+l+s=n, r,s>0 

with the boundary condition Gq = 1. When the product is not specified, it is taken to be 
the concatenation product *i . In this section, we call the generating function of moments 
Z(^) = X^n^"^"- This is because we will expand Z{^) in powers of q, and the coefficients 
ZkiO are not to be confused with the moments Gn, which are coefficients in an expansion in 
powers of ^ . Of course, q is a. bookkeeping device which is eventually set to 1 . 



3.4.1 Exact solution 

The loop equation for the Gaussian (|83|) may be solved since it is a quadratic equation 

a Z^{^)^ afZ^ - Z + 1 = 0. (84) 



The solution is 



z(() = '-y^ = Er.««^" (85) 

where r„ are the moments. Define Catalan numbers C„ by 

^ „ l-VT^ (2n)! 4" 

> CnX = With Cn = —r, — — TT ~ —F= as 77, ^ OO. (86) 

n\(n + l)\ ^J■n?^l 

Then the non- vanishing moments of the Gaussian 1-matrix model are 

(4a)" 

^2n = GnOi^ ~ —1= as n ^ oo. (87) 
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3.4.2 Approximate solution by deforming the product 



In ()83() . Dq is a derivation of sh = *q, not of cone = *i. So it is not a differential equation. But 
we can expand *i in a series in powers of q{= 1) around *o • Expanding Z{^) also in a power 
series in q , turns the loop equation into a sequence of differential equations in the shuffle algebra. 
At order , we get a nonlinear ODE for Zq{^) . Beyond that, we get a linear inhomogeneous 
ODE for Z].{^) in terms of In the end, q is set to 1. Let us illustrate this at 0{q^) 

and 0{q^). From section the expansion of *q around *o = sh is 

{F *, cm = (F *o cm - 1 e *o e *o {DoPm *o {DoGm + ■■■■ 



Moreover DqC = 1, so keeping only terms to 0{q), 

1 
2 



{Z *q (,)*qZ = (Z *o C - *0 DqZ *o ^ *0 DqO *q Z 



Z *o ^ *o Z - *o ^ *o Do{Z *o *o DqZ - -£,*oC *o DqZ *o Z 



Z *o Z - ^ 



2^ *o(,*o Z *o DqZ + C*o£.*o(,*o DqZ *o DqZ 



So the LE are 



DqZ = a 



C*oZ*oZ -^1^2^*0 C*oZ *o DqZ + ^ *o C *o C *o DqZ *o DqZ^ + 0{q^ 



(.90) 



Suppose Z{$^) = Zq{^) + qZi{^) + q^Z2{i) + . . .. Comparing coefficients of q^ and q^ we get 

— DqZq = Zq *o Zq 

a 

—DqZi = 2^ *o Zq *o Zi - - ^2^ *o ? *o Zq *o DqZq + C *o C *o *o DqZq *o DqZq^ . (91) 

So we have a non-linear ODE for Zo{S,) , and linear in-homogeneous ODEs for Zk,k > 1. The 
boundary condition Z{0) = 1 becomes Zo{0) = 1, Zfc(O) =0, k>l. 

Zeroth order 0{q^). Replace concatenation by shuffle product: The ODE for Zq can 

be linearized by passing to the shuffle reciprocal of Zq{(^) 

Y{C)*oZoiO = l DoZo = -Zo*oZo*oDoY. (92) 
Y satisfies the inhomogeneous linear ODE DqY = —a^ with boundary condition 5^(0) = 1. So 

Y{0 = 1-^^*1^ = l-ae- (93) 
Taking the shuffle reciprocal, we get (using ^^"'^ = n\ -^^^'^ = n! 

zoio = (1 - f e *o = 1 + f e *o e + (f )'e *o e *o e *o e + • • • = E ^('^^v- (94) 

n=0 

So the generating function at order q^ is 

^(0 = E U (2^)! + 0(^7). (95) 

n=0 ^ ^ ^ 
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And the non-vanishing moments in this approximation are G2n = y^j (2?^)! + C>{q)- These 
are compared with the exact moments in the table below. 



Moments 


exact 


0{q^) 


G2 
G4 
Ge 
Gs 

G2n,n 00 


a 
2a2 

5a^ 

(4a)" 


a 
6a^ 
90a3 

25200"^ 
(f)"(2n)! 



Due to the (2n)!, the 0{q^) moments numerically exceed the exact moments. We have a crude 
zeroth order answer with the potential for calculating corrections. Of course, the gaussian is a 
trivial model to solve. The value of our method lies in its applicability to multi-matrix models 
for which no method of solution exists. 



3.4.3 Approximate solution by deforming the left annihilation operator 

Next, we expand Dq around Di so that the loop equation (|55|) becomes a sequence of differential 
equations with respect to cone. This leads to a different approximation compared to section 
13.4.21 where we used the deformed product. Here, the expansion parameter is p = 1 — g, which 
is eventually set to 1. Recall that the (7 -deformed annihilation operator is 

D,F{x) = ^ D\F{i) = D^FiO - ^^DfFiO + ^eoMO + Oip'). (97) 

k=l 

If we expand Z{$,) in powers of p, Z(^) = J2n^n{(,)P^^ then 

00 s (_-\]n 

D,z{o = T.p'Il 7^^^r'^s-n(e)r (98) 

and 

00 s 

Z{0 *1 i *1 =T.P''T. *1 e *1 Zs-nii). (99) 

s=0 n=0 

Comparing coefficients of p, we get a nonlinear ODE for Zq{^) and a sequence of I''* -order 
linear ODEs for Zs{C) in terms of the lower order ones Zs„i(^), • • • : 

s (_-\\n s 

The first couple of orders are (all products are concatenation products) 



Z)iZo(e) = aZo(Oe^o(0 

) 

(101) 



Z)iZi(e) - ^D?Zo(Oe = 2aZo(e)eZi(e) 
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Zeroth order: At 0{p^) we have to solve the ODE DiZq{C) = a Zo(0 *i C *i ^o(0 with 
(0) = 1 • The solution is the cone reciprocal 

t2 ^,2 t4 „,3 



1 



, ae a^e' ^ 

1 + + — ^ + — ^ + ^ + + • • • = E 



n=0 



^^".(102) 



l-iae*iC ' ' 2 ' 4 ' 8 ' 16 
The non-vanishing moments are thus 

G2n=[^y + 0{p). (103) 

These are compared with exact moments = C'no" ~ , i , in Table IllUI We see that at 
leading order, deforming the annihilation operator underestimates the moments. 

Next to lowest order 0{p^): At the next order in p = 1 — g we have an inhomogeneous linear 
first order ODE for Zi{(^) 

DiZiii) - ]^DlZom = 2aZo{i)iZi{i) (104) 
with boundary condition ^i(O) = 0. Now Y' + PY + (5 = has solution 

y(0 = -I'^ii) f Q{r])I{r])dr] where !{£,) = exp f P{r]). (105) 
Jo Jo 

y = Zi(0; P = -2aCZo{0; Q = -^^Zl^O; Zoiv) = /(e) = (1 - ^a^' • Thus, 

v5 tlO 



3ae^ + 81og(l-^<^) _ 25 g6 41 ^4 ^8 99^5 



4(1 - la^)^ 



^ n>2 ^ ^ ^ r=0 



r + 1 



n — r 



(106) 



To get the asymptotic behavior of moments for large n, let Zi{S,) = J2 G2n^ 

n~2 



2n 



a 



2 — (-^2n 



2n 



r + 1, 



r=0 



, n > 2 



7 

2n log n - (- - 27)n + 2 log n + 0{n^] 



,n^oo. (107) 



Recall that Z{C) = Zo{^) +pZi{^) + -- ■ and Zo{0 = En(f )"^^"- Combining, at 0{p) we have 
(after setting p = I) 

n-2 „ , 1 n 

+ 0(p^), n > 2 



G2 = 

G2n ~ 



3a 



+ 0{p-); G2„ 



r=0 



n — r 



7 

2n log n - (- - 27)n + 2 log n + ©(n"' 



+ 0{p^),n 00. 



This is to be compared with the exact moments 

(2n)! ^„ (4a)' 



■ 2?i 



a ~ 



n!(n + l)! 



n — > 00. 



(108) 



(109) 
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Going to the next to leading order in p has improved the agreement with the exact correlations. 
For large n, the next to leading corrections to G2n are bigger in magnitude than the 0*^ order 
G2n ■ The accompanying table summarizes the approximate correlations obtained by expanding 
the left annihilation around the full annihilation operator in powers of p = 1 — q. 



Moments 


exact 


0{p'>) 


0{p) 


G2 


a 


0.5a 


0.75a 


G4 


2a2 


0.25a2 


0.75a2 


Ge 


5a^ 


0.125a3 


0.6460^ 


Gs 


Ua^ 


0.0625a^ 


0.490a^ 


G2n,n 00 


(4q)" 
V Trn^ 


(f)" 


(ff(2nlogn) 



(110) 



3.5 Non-Gaussian 1-matrix models 

Recall that the 1-matrix loop equation H12|) for a polynomial action tr S{A) = tr J2h=i ^i^^ 
with Sm 7^ determines higher rank correlations Gm-i,Gm,Gm+i, - ■ ■ in terms of the lower 
rank ones Gq = 1, Gi, G2, ■ ■ ■ Gm-2 ■ Suppose we apply our approximation method here. At O*'^ 
order we replace cone by sh. Since left annihilation D^"^ = is a derivation of sh, the loop 
equation becomes a quadratically non-linear ODE in the commutative shuffle algebra 

m 

J2 ISiD'-'GiO = G(0 o ^ o GiC). (Ill) 
1=1 

However, for m > 2 (i.e. non-Gaussian models), the differential operator X]l=i lSiD^~^ is not a 
derivation of sh and our trick of passing to the shuffle reciprocal does not linearize this ODE. 
It can still be thought of as a set of recursion relations (use *o = 

m I 

Y,lSiGr+i-i= for r = 0,1,2,... (112) 

1 = 1 s + t + l=r ^' 

S,t>0 

which determine Gm-i, Gm, Gm+i, • • • in terms of Gi, G2, ■ ■ ■ Gm-2 '■ 

r = : SiGo + 2S2G1 + ■■■ + mSmGrn~l = 

r = l: SiGi + S2G2 + --- + SmGm = l. e.t.c. (113) 

Our approach does not lead to a significant simplification for non-Gaussian 1-matrix models. 
However, we observe that the passage to the limit g = (replacement of cone by sh) did not 
change the dimension of the space of solutions to the original loop equations. 

4 Approximation method for multi-matrix models 

Recall the multi-matrix LE (|30|) for the generating series of gluon correlations 5*G(^) = G(^).^*G(^) 
where 5* = X]n>o("' + l)'^'''^ "'"*^jn " " " ^ji ■ Products on the RHS are cone products, but Dj 
are not derivations of cone. So the LE are not differential equations. By analogy with 1-matrix 
models, two ways around this mismatch come to mind. We could p-expand Di around full 
annihilation Dj, which is a derivation of cone. Or, we could g-expand cone around sh, with 
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respect to which Di is a derivation. Both these turn LE into quadratically non-hnear PDEs at 
Qth Qj-(jgj. in p or In the former approach these are PDEs on the non-commutative concate- 
nation algebra, while in the latter case, they are PDEs on the commutative shuffle algebra. We 
focus on the second approach in section H?T] due to its similarity with deformation quantization, 
and briefly consider the first approach in section 14.21 Beginnings of a formalism to go beyond 
zeroth order are in section 1^31 

4.1 Multi- matrix LE at 0(g°) and the shuffle reciprocal 

At O*'^ order in we replace cone by sh. Then the factorized LE H3U|) become^^ 

5^G(e) = G(e)oeoG(o (114) 

with the boundary condition G{0) = 1. (|114|) is a quadratically non-linear PDE on the shuffle 
algebra. In general, the order of the PDE is one less than the degree of S{A) . If 5* is a 
derivation of sh, we can change variables so that H114|) becomes a linear PDE for the shuffle 
reciprocal of G(^) denoted F{(^), -F(^) o G(^) = 1. The shuffle reciprocal exists as a formal 
series since the constant term Go = 1 does not vanish. Moreover, since Gj are cyclic and shuffle 
product preserves cyclicity, Fj are also cyclic. Assuming 5* is a derivation of sh, 

S'{F{C)oG{^))) = =^ F oS'G = -S'F oG S'G = -GoS'FoG. (115) 

Putting this in 1)1141) we get GoS^F oG = Go^*oG. Shuffle multiplying by F o F reduces the 
LE to a system of inhomogeneous linear PDEs in the shuffle algebra 

S''F{i) = -C- (116) 

We call these shujfle reciprocal LE. We seek cyclically symmetric solutions to them. The LHS 
of ()116() is the same as in the LE (|3U|) with G replaced by its reciprocal F . The — sign due to 
inversion has been written on the RHS. The RHS, however is much simpler than in (|30)) since the 
quadratic factor in G(^) has been eliminated. For the zero-momentum Gaussian, Chern-Simons 
and Yang-Mills matrix models we get (from section [27f|) 

Gaussian C'^ DjF{i) = - f 
Ghern- Simons in e'^''[Dk, Dj]F{C) = -f 

Yang -Mills -- [^fe, A]]i^(0 = " (117) 
a 

Thus, we have used the derivation properties of these theories to effectively linearize the LE at 
order . We still have to solve these linear PDEs on the oo -dimensional vector space spanned 
by First we find a formula to recover G(^) from its shuffle reciprocal F[^). 

(FoG)(0 = l^ E I'jGK = di. (118) 

I=JUK 

We can solve these equations starting from Go = i^o = 1 • The first few equations are 

Fi + Gi = 0, Fij + FiGj + FjGi + Gij = 0, 
^'^Here, we use o for *o = sh to avoid subscripts. 
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Fijk + FijGk + FikGj + FjkGi + FiGjk + FjGik + FkG-ij + Gijk = 0, ... (119) 

Since each successive equation involves the next higher rank Gj only linearly, we need only solve 
a linear equation at each step. Thus for |/| > 0, 

Gi = - Yl PjGk (120) 

I=JUK,K^I 

expresses higher rank Gi in terms of lower rank ones and the reciprocal F . Iterating, 

\l\ 

Gi = Y,{-ir E Fj,Fj,...Fj„ for / / 0. (121) 

n = l I — IiUl2U---UIn 

Ik 7^1>y k 

/ = /i U /2 U • • • U 4^ Ii, - ■ ■ ,In are complementary order-preserving subwords of I. For 
example, Gi = —Fi, Gij = —Fij + 2FiFj and 

Gijk = —Fijk + '^{FiFjk + FjFik + FkFij) — GFiFjFk 

Gijki = —Fijki + 2{FiFjki + FjFiki + FkFiji + FiFijk + FijFki + FikFji + FuFjk) 

-Q{FiFjkFi + FjFikFiFiFjiFk + FjFuFk + FkF.jFi + FiFkiFj) + 2AF^FjFkFll22) 

This formula shows that the mapping to shuffle reciprocal (for series with non- vanishing constant 
term) is one-to-one. We don't lose any information in going from G{S,) to F{^) and back. Once 
we solve ()116() . for F{S^) we may straightforwardly recover Gj using H121|) . 

4.1.1 Solution of Gaussian multi- matrix model at zeroth order in q 

Consider the Gaussian multi-matrix model tr S{A) = ^ tr C^^AiAj with symmetric covariance 
C«i = C^^ . At 0**^ order in q , the shuffle reciprocal LE (|117jl are 

DkFiO = -Ckji\ (123) 
where Ckj = Gjk is the matrix inverse of G^^ . We seek a solution of (|123|1 of the general form 

F{0 = 1 + Fi,e + Fi,i,e''' + ■■■+ i^n-i„f + • • • , (124) 

where Fj are cyclically symmetric. Gq = 1 fixes Fq = 1. Substituting in ()123() using = 

^h^i2-in get 

Fi + Fii.e' + Fii,i,e''' + ■■■+ Fii,...i,^e'-'" + ■■■ = c,,e- (125) 

Comparing coefficients of words we read off the solution 

Fi = 0, Fij = -Gij, Fi^...i^ =0 for n > 3. (126) 

The solution is a quadratic polynomial -F(^) = 1 — Gij^"^^ . Using (|121|) we get 

Gq = 1, Gi = 0, Gij = —Fij = Gij, Gijk = 0, Gijki = 2{GijGki + GikGji + GuGjk}, ■ ■ ■ (127) 

Thus, for the Gaussian multi-matrix model, the linear equations (|123|) for shuffie reciprocal, 
along with the boundary condition Fq = 1 have a unique solution. Comparing with exact 
moments from the planar Wick theorem, Tq = 1, Tij = Gij, Tijki = GijGki + GuGjk, ■ ■ ■ , we 
see that the approximation is an over estimate (as we found in the 1 -matrix example in section 
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4.1.2 Chern-Simons matrix model at zeroth order in q 

Consider the zero-momentum limit of 3d Chern-Simons(CS) gauge theory. This corresponds to 
the 3 -matrix model with action tr S{A) = 2iK tr Ai[^2)^3]- Such an action also results from 
considering terms in that are linear in momentum. The O*'^ order CS loop equation (|117() 
for the shuffle reciprocal F{(^) is 

if^e'^''[Dk,Dj]F{0 = -e or 2iK e'^^DkDjFi^) = . (128) 

We seek a solution to (|128() among formal series F{^) = Fj^^ with cyclic coefficients Fj satisfying 
Fj = Fj . Ean. (|128|) is an inhomogeneous 2"^^^ order linear PDE in an infinite dimensional space 
spanned by the words Fq = 1 does not suffice to fix a solution. For example, Fi are 
undetermined, since ^* is annihilated by the LHS. Inserting Fi^^ into ()128|1 gives 

2iK e''^'^Fi,...i„e^-'" = -e. (129) 

The PDEs become linear equations for the coefficients Fj with \I\ > 2, 

n = 2 : e'^'^Fjk = 0; n = 3 : 2iK e'^^Fjki = -5\; and n > 3 : e^^^^^i?.^ ..^ = q. (130) 

Being a system of inhomogeneous linear equations, the general solution is the sum of a partic- 
ular solution and the general solution of the corresponding homogeneous system. A particular 
solution with minimal number of non- vanishing Fi is 

Fo = 1, Fi23 = i^23i = i^3i2 = i^3*2i = ^2*13 = ^132 = ^ and = V other /. (131) 

To see this we need only consider n = 3, where despite appearances, after accounting for cyclic 
symmetry, there are only a pair of independent equations, the real and imaginary parts of 

i^321 -i^l23 = 7^. (132) 

By hermiticity, F^2i = -^123 or 3^-^321 = 3^-^123 and 9-F321 = — 9i^i23- Since k, is real, the real 
part of the above equation is an identity, so the real part 5R-F123 = 3^-^321 is left undetermined, 
and we can set it to zero. Its imaginary part gives 9-Fi23 = ^ , which is the advertised particular 
solution. For this particular solution the gluon green functions at order can be non-trivial 
only if their rank is divisible by 3 . For example. Go = 1 , 

Gi = -Fi = 0, Gij = 0, G123 = G231 = G312 = G321 = Gi32 = G213 = — , 
Gijk = otherwise , G113322 = -8G132-F132 = -77-2", G112233 = 0, etc. (133) 

Let us now consider the general solution to the inhomogeneous linear equations H13U|) . It is 
straightforward to see that they have infinitely many solutions, since the corresponding homo- 
geneous equations e**^*2Fj^j2...j„ =0, n>2 do. Indeed, any tensor Fi^i^i^...i^ that is symmetric 
under interchange of a pair of adjacent indices is a solution to the homogeneous equations. By 
cyclic symmetry, the two indices can be chosen as ii and ^2 • Then such an Fi^i^i^...i^ is annihi- 
lated due to antisymmetry of e**i*2 _ Even after imposing hermiticity and cyclic symmetry, this 
will leave an infinite number of homogeneous solutions, for example any totally symmetric real 
tensor Fi is automatically cyclically symmetric and satisfies Fj = Fj. To get an idea of how 
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many solutions there are among tensors of a fixed rank, consider eacli rank individually since 
the equations do not mix tensors of different rank. For n = 1, we do not have any LE, but 
hermiticity implies that Fi,F2,F^ are three arbitrary real quantities. For n = 2, e^^^Fjk = 
does not impose any condition on -Fii, ^22, ^33 ) which are real by hermiticity, and says that 
-^12 5-^23 and -F31 are symmetric tensors, which must again be real. For n = 3, as we saw earlier, 
2iK e^^^Fjki = —6\ is just the single condition 9-Fi23 = After accounting for cyclicity and 
hermiticity, there are 11 independent components of Fij^. The 10 undetermined components 
can be taken as the real numbers 

-f222, -^333 ) 3^-^123 5 -^"122, -^233) -P3II1 -^133) -^211, -P322- (134) 

For n = 4 , accounting for cyclic Fj , there are only 9 conditions 

-^"1123 = -^1132 = -^"1213 5 -^2231 = -P2213 = -^2321, -^3312 = -^3321 = -^3132, 

-^"1212 = -^1122, -f2323 = -^2233; -^3131 = -?^3311- (135) 

But there are c{n = 4, A = 3) = 24 independent cyclic symmetric fourth rank tensors (see 
appendix^. Thus we have a large space of homogeneous solutions among fourth rank tensors. 
A similar situation continues for n > 4. The LE at order (|128j) . though linear and easy to 
solve, have infinitely many solutions. As explained in section this is true of the original LE 
and is not an artifact of our approximation scheme. It remains to see if the additional equations 
obtained in 12.31 fix this shortcoming. 

4.1.3 Yang-Mills multi-matrix model at zeroth order in q 

Consider the Yang-Mills matrix model with action tr S{A) = - ^ tr [^j, ^j] [At, Ai\g^^ g^^ . The 
LE for the shuffle reciprocal F{^) of the moment generating series G(^) at zeroth order in q are 

g"'g^\Dj,[Dk,Di]]F{i) = ae for 1 < i < A. (136) 

Interesting special cases are A = 4, 2 which correspond to the zero momentum limit of 4 and 2 
dimensional large- Yang-Mills theory. For A = 2 and a flat Euclidean metric g"^^ = 5^^ , the 
matrix model action is tr S{A) = — ^ tr [^1,^2]^- l|136j) are a system of A inhomogeneous 
3rd Qj-fjgj- linear PDEs for F{^) = F^^j which is normalized to = 1- Fi must be cyclically 
symmetric. We need additional conditions to fix a solution since any quadratic polynomial 
is annihilated by the LHS, so Fi and Fij are not fixed by H136|) . Let us assume the metric 
gtj — ^tj g^j^j j^q|- niake a distinction between lower and upper indices, with repeated indices 
being summed. Then p36|) becomes (using the short-hand Dijk = DiDjDk) 

(2Djij — F)jji — Dijj)F{^) = a =^ (2Fjiji^...i^ — Fijji^...^^ — Fjjii^...i^)^ * " = a i^j(137) 

Comparing coefficients we get these conditions 

n = 4 2Fjijk - Fijjk - Fjjik = a 6ik V i,k 

^ ^ 5 =^ '^■^jijii---in -^ijji4---in ^jjii4---in — ^ ^ ^) ^4 ■ ■ ■ ^n- (138) 

The condition for n = 3 is an identity for cyclically symmetric tensors, so we drop it. These are 
infinitely many linear equations for the tensors Fj . A major simplification is that the equations 



27 



do not mix tensors of different ranks, i.e. the matrix defining the system is block diagonal with 
all blocks finite dimensional. Let us specialize to the simplest non-trivial case of the A = 2 
matrix model. We will show that a particular (cyclically symmetric) solution is 

ex 

Fq = 1, F1122 = -^2112 = -?^22ii = -^1221 = --^ and the remaining Fj = 0. (f39) 

The only non-trivial part of this particular solution involves the rank n = 4 tensors. The 
equations for the rest are homogeneous and they can be set to zero. For n = 4 we need to find 
a solution to 2Fjijk — Fijjk — Fjjik = a 6ik- These look like four equations, 

2F2121 — -^1221 — -?^22ii = a, 2F2122 — -^1222 — -P2212 = 0, 

2i^l212 — -^2112 — -?^1122 = a, 2F1211 — -F2III — -?^1121 = 0. (140) 
But there is only one independent non-trivial condition after accounting for cyclic symmetry 

F1122 — i^l212 = (141) 

Thus we see that Fq = 1,-Fii22 and cyclic permutations = —a/2 and all other F/ = is a 
particular solution. The gluon green functions at order are obtained via the shuffle reciprocal 
(|121|) which imply that non- vanishing correlations have rank divisible by 4, for example, 

Go = 1, Gi = Gij = Gijk = 0, G1122 and cyclic = — , and other Gijki = 0, etc. (142) 

Now comes the harder question of the general solution of the homogeneous linear system 

n — 4 =^ IFjiji^ ^ijjk -^jjik — ^ ^) ^ 

n ^ 5 =^ 1Fjiji^...i^ Fijji^...i^ Fjjii^...i^ — V 1,14. ■ ■ ■ in (143) 

For n = 4, as we saw before, there is only one non-trivial equation -F1122 = -^1212 • But there are 
c(n = 4, A = 2) = 6 independent cyclically symmetric rank 4 tensors (see appendix ^ which 
can be taken as F2222, -^"1222, -^1122, -^"1212, -^1112, Fim . Hermiticity F/ = Fj implies that 
all of them are real since reversal of order of indices can be achieved by cyclic permutations in 
each case. Thus the general solution for rank four tensors assigns 5 arbitrary real parameters to 
F2222, F1222, F1122 = F1212, F1112 and Full. 

For n = 5 , once the dust settles, there are only two non-trivial equations 

F11122 = F11212 and F11222 = F12122 (144) 

Taking account of cyclic symmetry, there are c(5, 2) = 8 independent rank 5 tensors, which 
can be taken as F22222, F12222, F11222, F12122, F11122, F11212, F11112 and Fum . In general these 
are complex, but hermiticity and cyclicity imply they are all real. Thus we have two linear 
constraints on 8 real parameters and therefore a six real-dimensional space of solutions to the 
shuffle reciprocal LE for rank 5 tensors: 

F22222, F12222, F11222 = F12122, F11122 = F11212, F11112, Finn (145) 

are freely specifiable real quantities. 

^^Recall that n = 3 was identically satisfied by cyclically symmetric tensors. 
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This abundance of solutions continues to hold for n > 6. It is easy to see that the homoge- 
neous linear equations "^Fjiji — Fijjj — Fjjij = have an infinite number of solutions. Observe 
that any tensor that is totally symmetric in any three adjacent indices^^ satisfies this equation. 
In particular, totally symmetric tensors are an infinite class of solutions. The underdetermined 
nature of the linear equations for F((^) is not an artifact of our approximation scheme. It is 
already true of the full LE as shown in section 12.21 It remains to implement the additional 
conditions (j25]) to see if they select a solution. 



4.2 LE with deformed left annihilator and the concatenation reciprocal 

We also have the option of approximating the LE (|3fl|) by replacing left annihilation Di by full 
annihilation Dj at zeroth order in an expansion in powers of p = 1 — q. Since Dj is a derivation 
of cone, this again turns the LE into non-linear PDEs, but this time on the non-commutative 
free algebra. As before, it is possible to convert the non-linear PDEs into linear PDEs by passage 
to the concatenation reciprocal. Recall that 5* = X^nl'^ + ^)^''^ ''"^F)j^ ■ ■ ■ Dj^^ . When Dj is 
replaced by Dj , we denote the resulting differential operator 

S* = ^(n + 1)5^--^"*D,„---D,,. (146) 

n 

Moreover, assume couplings are such that S* is a linear combination of iterated commutators 
of Dj and therefore a derivation of cone. This is the case for the Gaussian, CS and YM matrix 
models or any linear combination thereof. At zeroth order in p, the LE become 

S*G(0 = GiOeCiO- (147) 

Now we'd like to use the same trick as before and turn this into a linear equation for the cone 
reciprocal of G(^) . Though cone is non-commutative, left and right concatenation reciprocals 
of G(^) are both equal. Let GR = 1 and LG = 1. Multiplying the first equation by L from the 
left and using the second, we get R = L. So let F{^) be the unique two-sided cone reciprocal^'^ 
of Assuming is a derivation of cone, S'{FG) = {S'F)G + F{S'G) = 0. This turns 

(|147)) into a linear equation for the cone reciprocal 

S'F{() = -f . (148) 

Inserting F(^) = F/^^ into (tTlHl) gives linear equations for coefficients Fj . Once Fj are deter- 
mined, we recover Gj at zeroth order in p using the following formula for cone reciprocal. 

{FG)i = 6j^Q ^ Go = l and ^ji^^F/^G/^ = Gj = - ^ F/.G/^ for |/| > 0. (149) 

Iterating this, we solve for Gj 

\i\ 

Gj = J2{-ir E FhFh---Fi,, for 1/0. (150) 

n=l I=Ill2---I-n 

7fc ^ Vfe 



^^By cyclic symmetry, those three indices can be taken as the first three. 
20^0 = Go = 1 
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For example, the first few gluon correlations are 



Go — 1; Gi — —Fi] Gij — —Fij + FiFj] Gijk — —Fij^ + FijFf^ + FiFj^ — FiFjFj.; . . . (151) 



Thus cone reciprocal is a 1 — 1 map. However, unlike shuffle reciprocal, it does not preserve 
cyclicity. Though we do not lose any information in the passage from G to F , the cyclic property 
of Gj gets slightly garbled when expressed in terms of Fj . For example, cyclic symmetry of 
Gijk implies the relation FijF^ — Fijk = FjF^i — Fjki- Thus, we should look for solutions to 
H148() among Fj that lead to cyclically symmetric Gj 's. This makes identifying the appropriate 
solutions of (|148j) potentially harder than for the corresponding shuffle reciprocal LE (|116j) . 
There is another reason why the concatenation reciprocal LE (|148|) are a potentially harder 
infinite linear system to solve than their shuffle reciprocal counter part H114|) . Left annihilation 
acting on a monomial produces a monomial [DjF]i = Fji . But due to its democratic nature, full 
annihilation produces a linear combination of monomials [DjF]/ = Sj^^^Fj^^jj^ . Thus the matrix 
defining the system of linear equations for Fj would be less sparse than before. Nevertheless, 
the moral is that replacing Dj by Dj at O*'^ order in an expansion in p allows for an effective 
linearization of the LE provided the action has the derivation property. 

4.3 Formalism for multi-matrix models beyond zeroth order 

At 0{q^) our approximation amounted to replacement of non-commutative cone by commuta- 
tive sh in the LE. This is like approximating the associative product of operators in quantum 
mechanics by a commutative product of functions on phase space. To go beyond this, we need 
a formula expressing eonc as a series around sh, hy analogy with the Moyal * -product formula 



for the symbols of operators (here Weyl ordered) in quantum mechanics. The first non-trivial 
term in such a formula involves the classical Poisson bracket. So one strategy is to look for a 
natural Poisson bracket on the shuffle algebra. However, there are differences from the usual 
situation where Heisenberg equations are approximated by Hamilton's equations. While the 
Heisenberg equations of quantum mechanics involve commutators of the associative product, 
the LE directly involve the associative concatenation product and not its commutator. Another 
difference from the usual situation in deformation quantization is that we know the product 
at both q = and q = 1 whereas one usually knows the product only at h = 0. Once we 
have such a formula, then as we did for 1 -matrix models (section l3.4.2|l . we would expand the 
generating series of gluon correlations G{() = EkLoG^'^HOQ^ 

in a power series in q and find 

equations for the G*-*^^ (^) order by order in q , starting from the 0*^ order equations for G^^^ (^) 
of section ITTl However, the situation for multi- matrix models is substantially more complicated 
than for the 1 -matrix models of section |31 This is because cone is non-commutative while it 
was commutative in the single-matrix case. 




where 



etc (152) 
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4.3.1 g-Deformed product and Poisson bracket on shuffle algebra 

We exhibit a 1 -parameter family of associative products *q that interpolate between commu- 
tative shuffle *o and concatenation *i . It reduces to the g -product for a single generator 
introduced in ((70)) and is defined as {F *q G){S,) = [F *q G]i(,^ where^^ 

[F*qG]i= Yl p''^''^'''^FjGk and p = 1 - q. (153) 

JUK=I 

The (two-word) crossing number xi^'i K) of the ordered triple {I; J, is the minimum 
number of transpositions of elements of J and K in order to transform JK into / when J and 
K are order-preserving sub- words of /. For example, 

X{ijk;i,jk) =0, x{ijk;ik, j) = I, x{ijk; jk,i) = 2. (154) 

For q = 1 {p = 0) , this formula reduces to cone. For, the only term that contributes is the one 
with xi^i JjK) = i.e. no crossings, so I = JK . Then 

{F*iG)i = 6i''FjGK. (155) 

If q = {p = 1) , then p^i^'-^'^) = 1 independent of the crossing number and all terms contribute 
equally giving back shuffle 

{F*oG)i= J2 ^J'^K. (156) 

I=JUK 

Examples: For g 7^ 1, *g is non-commutative in general. The first few terms in the g-product 
of a pair of tensors are {F *g G)o = FqGq , 

{F*gG)i = FiGo + FoGi, {F *g G)ij = F^Gij + F.Gj + pFjGi + Fi.Go, 
{F *q G)ijk = FoGijk + FiGjk + pFjGik + p'^FkGij + FijGk + pFikGj + p'^FjkGt + FijkGo 
(F *g G)ijki = FoGijki + {FiGjki + pFjGiki + P^FkGiji + p^FiGijk) 

+{FijGki+pFikGji+p^FuGjk +p^FjkGu+p^FjiGik+p'^FkiGij) 
+{FijkGi + pFijiGk + p^FikiGj + p^FjkiGi) + FijkiGa. (157) 

Associativity: We show that the g-product is associative 

{{F*gG)*gH)i = {F*g{G*gH))i= ^ ^'^'^^ Fj i^L • (158) 

I=JUKUL 

We first checked explicitly that associativity holds for \I\ < 3 by writing out all the terms, but 
it was very tedious to go further. Instead, we write 

{{F*G)*H)i = p''^''''''\f*G)jGk= E p'^^''''^*''''V'^^''^*'''"*'^i"LG'M^i^ 

I=JUK I=LUMUK 
I=JUKUL 

{F*{G*H))j = £ pX{I;J,KUL)+x(KUL;K,L)p^^^jj^ (^59) 

I=JUKUL 

^^To avoid too much clutter we will occasionally drop the subscript in and indicate it by *. 
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where / = JUKUL is the condition that J, K, L are complementary order-preserving sub-words 
of / . Since F, G, H are arbitrary and so is p , associativity requires the equahty of the sums of 
crossing numbers 

x{I;JUK,L)+x{JUK;J,K) and x{I; J, K U L) + x{K U L; K, L) (160) 

for each / and any (fixed) choices of J,K,L,J\JK and KUL satisfying I = JUKUL. In fact, 
these two sums of (two-word) crossing numbers are equal to the (three-word) crossing number 
xil] J, K, L) that has a simple meaning. x(I; J, K, L) is the smallest number of transpositions 
needed to transform JKL into / where J, L are order-preserving sub- words of / . For 
example suppose / = abed, J = d,K = c,L = ab,J U K = cd and K U L = abc. Then 
x{0'bcd;cd,ab) + x{cd]d,c) = A + 1 = b while xi'Acd^d, abc) +x{ohc;c,ab) = 3-1-2 = 5. Similarly, 
if / = abed, J = b,K = ad,L = c, J U K = abd and K U L = aed. Then xi^^bcd; abd, c) + 
x{abd;b,ad) = = 2 while x{(^bed;b,acd)+x{acd;ad,e) = 1-1-1 = 2. Thus, associativity just 
says that there are two different ways of calculating the three- word crossing number xi^'i L) 
when / = J U K U L . This gives the simple formula 1)1581) for the *^ product of three series, 
which makes associativity manifest. 



Reduction to one generator: When we reduce to a single generator in the above examples 
(|157)) . the formulae agree with those obtained earlier ()77|) using the Gauss binomials. More 
generally, we can see from the definition of the Gauss binomials H71|) that 

q I^JUK 
\J\=r 

Thus, the above formula for the g-product reduces to the one for a single generator. 




Poisson Bracket: It may help to find a Poisson bracket on the shuffle algebra that serves as a 
first approximation to the g -commutator. The (7 -commutator is 

{[F,G]q)i^{F*qG-G*qF)j= J2 {l-qf^''''''HFjGK-GjFK). (162) 

I=JUK 

For small q, -l{[F,G]g)i = Ei=juK xil; -l K){FjGk - GjFk) + 0{q). So let us define the 
bracket {F,G} = {F,G}iC^ by 

{F,G}i = -hm-{[F,G]g)i= J2 x{I;J,K){FjGk-GjFk). (163) 

^ I=JUK 

It is clearly bilinear and anti-symmetric. The first few examples with lowest |/| are 

{F,G}o = 0; = 0; {F,G}^J = FjG., - GjFf, 

{F, G}ijk = FjGik + 2FkGij + Fi^Gj + 2FjkGi — (F ^ G); 
{F, G}ijki = FjGiki + 2FkGiji + SFiGijk + FikGji + IFuGjk + 2FjkGii 

+3FjiGik + 4.FkiGij + Fi^iGk + 2F,kiGj + SF^kiGi - {F ^ G). (164) 

It satisfies the Jacobi identity since the g-product was associative. 

{{F, G},H}+ {{H, F},G} + {{G, H},F}= 0. (165) 
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This can also be checked exphcitly. For example, the first non-trivial case is 

{{F, G}, H}i,k = 2{FiGjHk + FkGjHi - F^GkH., - F^G^Ht). (166) 

Upon adding its cyclic permutations, the Jacobi identity is satisfied. Moreover, the Leibnitz rule 
(with respect to sh = o = *o) 

{FoG,H} = Fo{G,H} + {F,H}oG (167) 

is also satisfied due to the corresponding identity for the g-commutator. Thus {...} is a Poisson 
bracket on the commutative shuffle algebra. 

In order to be practically useful in going beyond the 0**^ order solution of the LE, we need 
a g-expansion for *g around *o = sh involving left annihilation Dj. For small q, 

iF*gG)i = J2 il-qV^''''''^FjGK= E I'jGK-q E xil; J, K)FjGk + Oiq'') 

I=JUK I=JUK I=JUK 

= {F*oG)i-q E x{I;J,K)FjGK + Oiq^) 

I=.JUK 

^ iF*,G-F*oG), ^ x(/;J,i^)F.G^. (168) 

g-*0 —q ^ 



I=JUK 



For example, 



(F*^G-F*oGk 



lim^ "i = F.Gi 

9-0 -q ^ 



{F*gG_J^*oG)ijk ^ F,G,k + 2FkG,,+FikGj+2F,,G,. (169) 



Our aim is to express this 0{q) contribution to F *g G in terms of Df and *o . But we are yet 
to find such a formula that generalizes H79|) and hope further investigation will reveal it. 



4.3.2 g-Deformed annihilation 

There is one parameter family of annihilation operators Dj that interpolates between left 
annihilation Dj and full annihilation Z)j . For a single generator, it was defined in H81|) as 
{DgG)n = (1 + g + (?2 + . . . + q-)Gn+i . By analogy we define [D]G]j = 6'j''^q\'^\Gj,jj, , i.e. 

[F^jG]ii---i„ = Gji^...i^ + qGij^ji2...i„ + q Gi-ii2ji3---i„ + •••+<? (170) 

We pick up one more power of q as the annihilation operator travels through each index of the 
tensor from left to right. It is easily seen that 

lim [L>?G]/ = Gji and lim [DfG]/ = dj^^^Gi^ji^ (171) 

q — s-0 q — ^1 

reproduce left and full annihilation which are derivations of sh and cone. To make the LE (|30|) 
differential equations with respect to cone, we want to expand Dj around Dj in powers of 
p = 1 — q and finally set p = 1. Recall that for 1 -generator (|82|) . 

DM^) = E ^i-PO'-'D',G{0 = D,G{0 - ^^DjGiO + ^enfGiO + 0{p').{172) 



33 



For several generators, 




+ 3G, 



Gji 



«l«2«3J*4---«n 



+ ■■■ + G^ 



+ ■■■ + 



n(n — 1) 



2 



P Gi 




H + nGi 



(173) 



Drawing inspiration from (|82)) we would like to recognize the coefficients of powers of p as 
combinations of full annihilation and some multiplication operator acting on G. However, we 
have not yet succeeded in this. 

5 Discussion 

Despite their formidable reputation, the loop equations(LE) of a large- multi-matrix model 
show much simplicity and structure when expressed in terms of gluon correlations Gj . Non- 
linearities are mild in the sense that in any equation, highest rank correlations appear linearly. 
So the LE are systems of inhomogeneous linear difference equations for correlations of a given 
rank with lower rank correlations appearing non-linearly as 'sources'. Solving these equations 
in the absence of additional structure would be tedious at best. But this is not possible because 
the LE are underdetermined in most interesting cases. We observed that there are additional 
equations involving the Gj that a naive passage from finite Schwinger-Dyson equations to 
large- A'^ LE misses. These equations have to do with changes of variables in matrix integrals 
that leave both action and measure invariant. However, we are yet to implement these additional 
constraints in detail to see whether they suffice to fix a unique solution to the LE. On the other 
hand, we saw that part of the difficulty in understanding the LE lies in the fact that they are 
not differential equations. Left annihilation does not satisfy the Leibnitz rule with respect to 
the concatenation product appearing in these equations. We proposed two schemes to remedy 
this situation by expanding either annihilation or product around one that is a derivation of the 
other. For the Gaussian, Chern-Simons and Yang-Mills models, it was possible to altogether 
eliminate the non-linearities of the LE and arrive at inhomogeneous linear PDEs at the zeroth 
order of these expansions. But the under-determinacy of the loop equations prevented us from 
picking a unique solution except in the case of the gaussian, where the two approximations 
were shown to give over and underestimates for correlations. This underscores the importance 
of better understanding the remaining constraints on Gj (section 12.3)) as well as any other 
conditions that would ameliorate the under-determinacy of the LE. In jl3j we hope to extend 
these algebraic and differential properties to matrix models with both gluon and ghost matrices, 
of the sort appearing in the gauge-fixed action of Yang-Mills theory. 
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A Cyclically symmetric tensors of rank n 



What is the dimension c(n, A) of the space of cychcally symmetric real tensors Gi^...i^ of rank 
(=number of indices) n if the indices can take the values 1 < ifc < A? The dimension of the 
space of all tensors of rank n is A". On the other hand, the space of symmetric rank n tensors, 
which is a subspace of cyclically symmetric tensors, is {^~^'^~^) dimensional. Thus 



'K + n-l 
n 



< c(n,A) < A" (174) 



For a A = 3 matrix model, |(n^ + 3n + 2) < c(n, 3) < 3". For a 2-matrix model, n + 1 < 
c(n, 2) < 2". The cyclic group of order n acts on rank n tensors Gi^...i^ by cyclically permuting 
indices. c(n. A) is the number of orbits. For example, if A = 2 and n = 4, the orbits are 

(G2222); (^1222 = ^2122 = ^2212 = ^2221); (^1122 = G2112 = G22II = G1221); 

(^1212 = G2121); ((^1112 = G2111 = G1211 = ^1121); (Giiii) (175) 

So c(n = 4, A = 2) = 6, significantly less than 2^. The cardinality of different orbits are not 
necessarily equal. Some other examples are 

c(n, 1) = 1; c(l,A)=A; c(2. A) = ^A(A + 1); c(3, 2) = 4; 
c(3, 3) = 11; c(4, 2) = 6; c(4, 3) = 24; c(5, 2) = 8. (176) 

It would be nice to have formula for c(n. A) , at least for the A = 2 matrix model. 



Note on hermiticity condition: Actually, the tensors G/ are complex numbers, so the 
real-dimension of the space of cyclically symmetric tensors of rank n is 2 c(n, A) . However, 
the hermiticity condition Gi^i^—i^ = G*^...^^^_^ halves this real-dimension to c(n. A). If reversal 
of indices can be achieved by a cyclic permutation (e.g. G1122 = G^22ii — ^2211 ) then the 
correlation is real. If I cannot be obtained from / via cyclic permutations, then hermiticity 
means that ^Gi = "RGj and = —QGj. For example 3^Gii23 = 5^6*3211 and 9Gii23 = 

—9^3211 . In either case, hermiticity halves the number of independent parameters in cyclically 
symmetric correlations of a given rank. 



B Concatenation, shuffle and their co-products 

By V , let us denote the infinite-dimensional complex vector space spanned by the monomial 
words in the A non-commuting sources ^* . A typical element is the formal series G(^) = 

G/^^ . V is the basic arena for our algebraic study of the loop equations^^ . 

The concatenation product cone : V (Si V ^ V denoted by juxtaposition, was defined in 
^^^"^ = 6j/ = ^^■^ . It has the structure constants c^jf = 6j/ . For A > 1 , cone is non- 
commutative. The vector space V , along with the concatenation product is the free associative 

^^A superior approach that makes cycUc symmetry of Gi manifest might be to consider the quotient by the 
relation ~ ^"^ if / is a cyclic permutation of J. Then a basis for V would consist of words where / labels 
orbits of the cyclic group action. In this paper we just allow all words and impose the condition that Gi be 
cyclically symmetric, by hand, so to speak. 
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algebra T on the generators ? , • " " i • It is the universal envelope of the free Lie algebra. The 
commutative shuffle product sh : V V ^ V was defined in (jlTj) . V , equipped with sh is the 
shuffle algebra. The shuffle product of monomials 

eoe = s'/e= E (177) 

IU,J=K 

leads to the shuffle structure constants s]^"' = \{I U J = K} \ . 

There is a natural inner product (., .) on V , for which form an orthonormal basis 

(^i,^J) = 5'^J or {Fje,GA') = FiGj5'^' = J2l'iGi. (178) 

I 

The Kronecker symbol 6^'"^ = 1 if / = J and otherwise. We can use the 'metric' 6^'"^ and 
its inverse 5jj to raise and lower indices. The inner product allows us to define co-products 
V V V . We call them co-concatenation A = s/i"^ and co-shuffle A' = conc^ . They are 
adjoints of sh and cone respectively. For three formal series F,G,H , we define A and A' by 

{F0G,A{H)) = {FoG,H) and {F (g) G, A' {H)) = {FG, H). (179) 

We define the structure constants of co-concatenation and co-shuffle as 

A(e'^) = sIm^ and A'(e^) = cIm^"^ ■ (180) 

We use the same letter c to denote the structure constants of cone and A' = cone^ because 
they are related by raising and lowering indices using the metric 6^''^ . The same goes for the 
letter s for the structure constants of sh and A = sh^ . The expressions for these are 

Cj,K = c^/'s''^6j,LdK,M = S'jK and s'j^j, = s^/' 6^'^ 5j,lSk,m = sf"" = \{I = J u (181) 
To obtain the co-shuffle structure constants Cj j^ , we use the definition of adjoint to get 

c^^P = 6'^'''6i,nSj,p = 6§p. (182) 
We use a similar procedure for the co-concatenation structure constants Sjj^ 

^ s£ d ' = Sl^m(> ^ =^ Sp,Q = ^L (> (>i^pd,iQ = Sj^^. (183) 

On formal series, co-shuffle A' = conc^ acts as 

A'F = [A'F]ij C^(S)^^ = Fui^ ® (184) 

In particular, A'(^-^) = ® and A'(^*) = 1 + 1 ^*). On formal series, co- 

concatenation A = sh'^ acts according to 

AF = [AF]j^Ki-^ (^i^ where [AF]j^k = E ^i- (l^^) 

I=JUK 

In particular, A(^^) = Y.i=juK C'^ <^ and A(f ) = ^» ® 1 + 1 
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C Bialgebra structures on V = Span{^^) 



V has two bialgebra (algebra + compatible coalgebra) structures. In one, the product {sh) 
is commutative while the co-product (adjoint of cone) is non-co-commutative. In the dual 
bialgebra, the product {cone) is non-commutative while the co-product (adjoint of sh) is co- 
commutative. 

To establish that shuffle and co-shuffle^'^ combine to define a bialgebra on F , we show that 
co-shuffle A' = cone^ is a homomorphism of the shuffle product 

A'(FoG) = A'(F)oA'(G). (186) 

Note that the LHS is 

M{FoG)= Yl FlGmA'{^'')= Y1 FlGmS,'®^-'. (187) 

LUM=K LUM=IJ 

While the RHS is 

A'(F)oA'(G) = wjA'(e^)oA'(eO = wXL^/,^(e''®e'')o(e'^^^^) 

= FKLGMNi^"" o^^')0{^^o^'')= Yl FklGmN^' ®^im 

KUM=I,LUN=J 

Comparing coefficients. A' = conc^ is a homomorphism of the shuffle product if 

E Pj^G-r^ = E Fl,l,Gm,m, V /i, h. (189) 

JlUj2=/l/2 LlUMl=Il,L2UM2=l2 

To prove this, observe that Ji may be uniquely decomposed as Ji = L1L2 with Li C Ii and 
L2 C I2 and similarly for J2, J2 = M1M2 with Mi C Ii and M2 C /2. Then we observe that 
every riffle-shuffle Ji U J2 = Iih arises from a unique pair of riffle-shuffles Li U Mi = 7i and 
L2 U M2 = I2 ■ This establishes that co-shuffle A is a homomorphism of sh . 

A similar argument shows that A = sh^ is a homomorphism of cone: A{FG) = A{F)A{G) . 

A(F)A(G)= E {AF)j,,j,{AG)j,,j,= ^ Fj.ui^Gi.uu- (190) 

J=hl3,K=l2l4 J=hh,K=l2U 

On the other hand, the LHS gives 

[A{FG)]j,K= E iFG)L= E Fl,Gl2= E ^/lU/^Gw^- (191) 

L=JUK LiL2=JUK J=Iil3,K=l2l4 

In the last equality, we used the unique decomposition J = IiI^,K = I2IA where /i,/2 C Li 
and /3,/4 C L2 as before. Thus we have shown that A = sh'^ is a homomorphism of cone. 

The unit element for cone isl, Fl = IF = F . The co-unit for co-concatenation is e : 
y — C. It picks out the constant term in a formal series e{Fj^^) = F^ = Fq. Just like 
co-concatenation, the co-unit is a homomorphism of cone 

e{FG) = {FG)o = FqGo = e{F)e{G). (192) 
^^This justifies the name co-shuffle for the adjoint of cone. 



37 



The unit element for shuffle too is 1 , (F o 1)/ = YjI=juk ^J^k — ■ The co-unit for co-shuffle 
is again e : y — > C . The co-unit e is a homomorphism of the shuffle product 

e{FoG) = {FoG)o= ^ FjGk = F^Gq = e{F) o e{G). (193) 

JUK=% 

To summarize, {cone, sh^ = A = co—conc, 1, e) defines a non-commutative but co-commutative 
bialgebra (algebra plus compatible co-algebra) structure onV = span{^^) . Similarly, {sh, conc^ = 
A' = co—sh,l,e) defines a commutative but non-co-commutative bialgebra structure on V. 
These two bialgebras are not independent. Structure constants of the product and co-product 
of one can be obtained from those of the other using the inner product 6^''^ on V . 

Remark: In addition to being a bialgebra T = (cone, A, 1, e) , is the universal envelope of the 
free Lie algebra. So it is a Lie algebra with the Lie product [^^, ^"^J = ^^"^ — ^"^^ . Does A define 
a Lie bialgebra (211 with respect to the commutator? No! On the one hand, A : V ^ V <^V is 
not skew-symmetric. Rather, its image lies within Sym(y (8) V) . 

A(e')= E ^''^^''= E e^C' = irA)ie), (194) 

I=JUK I=JUK 

where r(a<8)fe) = 6® a. Here we used the fact that if J and K are order preserving complementary 
subwords of /, then so are K and J. Furthermore, A is not a 1-cocycle for the free associative 
algebra. In order to be a 1-cocycle, it must satisfy 

A[F, G] = (adp 1 -Fl ® adF)A{G) - F ^G (195) 

for any F,G £ F . However, taking F = ^' and G = gives 

LHS = A[C,i^] + (196) 

and RHS = 2 x LHS ^ LHS . There may be some other skew-symmetric 1-cocycle A : T — > 
T ®T which defines a Lie bialgebra structure on the universal envelope of the free Lie algebra. 
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